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THE  INSTITUTE  FOR  APPLIED  TECHNOLOGY  provides  technical  services  to  promote 
the  use  of  available  technology  and  to  facilitate  technological  innovation  in  industry  and  Gov- 
ernment; cooperates  with  public  and  private  organizations  in  the  development  of  technological 
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plication of  radiation  to  the  solution  of  Bureau  mission  problems  and  the  problems  of  other  agen- 
cies and  institutions.  The  Center  consists  of  the  following  divisions: 

Reactor  Radiation — Linac  Radiation — Nuclear  Radiation- — Applied  Radiation. 

THE  CENTER  FOR  COMPUTER  SCIENCES  AND  TECHNOLOGY  conducts  research  and 
provides  technical  services  designed  to  aid  Government  agencies  in  the  selection,  acquisition, 
and  effective  use  of  automatic  data  processing  equipment;  and  serves  as  the  principal  focus 
for  the  development  of  Federal  standards  for  automatic  data  processing  equipment,  techniques, 
and  computer  languages.  The  Center  consists  of  the  following  offices  and  divisions: 

Information  Processing  Standards — Computer  Information  — Computer  Services  — Sys- 
tems Development — Information  Processing  Technology. 

THE  OFFICE  FOR  INFORMATION  PROGRAMS  promotes  optimum  -dissemination  and 
accessibility  of  scientific  information  generated  within  NBS  and  other  agencies  of  the  Federal 
government;  promotes  the  development  of  the  National  Standard  Reference  Data  System  and  a 
system  of  information  analysis  centers  dealing  with  the  broader  aspects  of  the  National  Measure- 
ment System,  and  provides  appropriate  services  to  ensure  that  the  NBS  staff  has  optimum  ac- 
cessibility to  the  scientific  information  of  the  world.  The  Office  consists  of  the  following 
organizational  units: 

Office  of  Standard  Reference  Data — Clearinghouse  for  Federal  Scientific  and  Technical 
Information  — Office  of  Technical  Information  and  Publications — Library — Office  of 
Public  Information— Office  of  International  Relations. 


1 Headquarters  and  Laboratories  at  Gaithersburg,  Maryland,  unless  otherwise  noted;  mailing  address  Washington,  DC.  20234. 
- Located  at  Boulder,  Colorado  80302. 

:1  Located  at  5285  Port  Royal  Road.  Springfield,  Virginia  22151. 


NATIONAL  BUREAU  OF  STANDARDS  REPORT 


NBS  PROJECT 

4314561 


June  1971 


NBS  REPORT 

10  432 


A SEARCH  AND  RESCUE  SIMULATION  MODEL  FOR  THE  UNITED  STATES  COAST  GUAR 

VOLUME  III 

PROGRAMMER  LEVEL  DOCUMENTATION  FOR 
"PREPROCESSOR” 


by 

S.  S.  Karp,  L.  K.  Cummings,  M.  D.  Maltese,  A.  L.  Weber 


Sponsored  by 
U.  S.  Coast  Guard 


NATIONAL  BUREAU  OF  S' 
for  use  within  the  Government, 
and  review.  For  this  reason,  tt 
whole  or  in  part,  is  not  authc 
Bureau  of  Standards,  Washing 
the  Report  has  been  specilicall; 


IMPORTANT  NOTICE 


Approved  for  public  release  by  the 
Director  of  the  National  Institute  of 
Standards  and  Technology  (NIST) 
on  October  9,  2015. 


sss  accounting  documents  intended 
i subjected  to  additional  evaluation 
; listing  of  this  Report,  either  in 
ib  Office  of  the  Director,  National 
by  the  Government  agency  for  which 
copies  for  its  own  use. 


<NBS> 

U.S.  DEPARTMENT  OF  COMMERCE 

NATIONAL  BUREAU  OF  STANDARDS 


PREFACE 


This  volume  is  one  of  a series  which  documents  a Search  and 
Rescue  Simulation  Model  for  the  United  States  Coast  Guard.  The 
material  reported  in  this  documentation  was  developed  by  an  inter- 
disciplinary team  at  the  National  Bureau  of  Standards  with  representa- 
tion from  the  U.S.  Coast  Guard  under  MIPR  Z-70099-0-01935. 

The  complete  documentation  is  comprised  of  the  following: 

Volume  I Executive  Level  Documentation 

Volume  II  Analyst  Level  Documentation 

Volume  III  Programmer  Level  Documentation  for  "PREPROCESSOR" 

Volume  IV  Programmer  Level  Documentation  for  "OPSIM" 

Volume  V Programmer  Level  Documentation  for  "POSTPROCESSOR" 

Appendix  A Flow  Charts  for  Programmer  Level  Documentation 
Appendix  B Program  Listings  for  Programmer  Level  Documentation 
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I . Introduction 

The  objective  of  the  Preprocessor,  PREPRO,  is  to  prepare  a set 
of  Search  and  Rescue  cases  to  be  simulated  in  OPSIM.  Figure  1 is  a 
macro  level  description  of  PREPRO.  Here  it  is  shown  that  the  SAR 
Assistance  Reports  go  through  a series  of  cleaning  routines  to  produce 
the  OPFAC  FILE*  which  contains  a set  of  cases  organized  by  OPFAC. 

Multi -unit  cases  are  assembled  via  a sorting  program  called  MUTAPE  to 
produce  a CASE  FILE.  To  include  the  C-130  interdistrict  cases,  a 
program  called  MUC130  merges  the  C-130  cases  with  the  CASE  FILE  to 
form  the  PCP  CASE  FILE.  MUTAPE  and  MUCL30  are  called  the  CASE 
ASSEMBLY  PROCESSOR. 

Once  the  PCP  CASE  FILE  is  prepared,  the  remaining  major  steps  of 
PREPRO  can  be  carried  out.  The  first  step  calculates  the  case  para- 
meters from  the  input  PCP  CASE  FILE;  the  software  designed  to  accomplish 
this  task  is  designated  the  Parameter  Calculation  Program  (PCP) . The 
PCP  produces  the  SARSIM  CASE  FILE  for  input  to  the  second  step,  the 
Demand  Generator,  DEMGEN,  or  the  HIST  program  which  orders  the  cases 
chronologically.  DEMGEN  or  HIST  produces  a CASE  TAPE*  the  exogeneous 
event  tape  input  to  OPSIM. 

Several  options  are  offered  to  the  user  in  preparing  the  final 
CASE  TAPE: 

*7TFTLEtt  refers  to  a collection  of  SAR  cases  not  necessarily  arragned  in 
chronological  order.  "TAPE"  refers  to  a FILE  in  which  the  cases  are 
arranged  chronologically. 


MAJOR  STEPS  IN  PREPRO 


DEMGEN 

GENERATE 

DEMAND 


PCP 

CALCULATE  CASE 

p*  1 

PARAMETERS 

S C ADC TM  1 .. 

"1  CASE  FILE  J 

HIST 

ORDER  CASES 
CHRONOLOGICALLY 


• The  user  may  use  the  historical  dates  and  times  of 

occurrence  of  all  cases  5 and  the  associated  case  para- 
meters with  the  HIST  option. 

© The  user  may  randomly  generate  the  dates  and  times  of 
occurrence  according  to  hourly  exponential  distributions , 
and  the  full  case  descriptions  relative  to  the  generated 
time  period. 

« The  user  may  increase  or  decrease  the  case  load,  and  choose 
to  do  this  relative  to  the  specific  case  type,  the  specific 
station,  or  across  all  cases  in  the  entire  district. 

© The  user  must  specify  the  scenario  limits -that  is,  the  exact 
time  span  over  which  the  simulation  will  be  run. 

© The  user  has  the  option  to  validate  his  generated  occurrence 
times  against  the  input  values  of  the  arrival  rates. 

Normally,  the  mode  is  to  run  the  CASE  ASSEMBLY  PROCESSOR,  and 
PCP  once,  generating  the  final  SARSIM  CASE  FILE  with  single  unit,  multi- 
unit, and  C-130  case  descriptions,  in  that  order.  Many  runs  can  be 
made  by  DEMGEN,  depending  on  the  user  selected  option. 

The  next  sections  discuss  the  CASE  ASSEMBLY  PROCESSOR;  PCP  and 
DEMGEN. 

II.  Case  Assembly  Processor 

The  case  Assembly  Processor  is  two  FORTRAN  V programs  which 
produce  the  tape  that  is  input  to  PCP. 

MUTAPE  assembles  the  multi-unit  cases  from  the  district  tape 
(OPFAC  FILE)  by  gathering  together  from  the  participating  OPFAC's 
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for  a single  case  all  of  the  C --cards  (SAR  Reporting  Command  data) 
and  associated  D-cards  (Assisting  Resource  data) . These  multi -unit 
cases  are  gathered  at  the  end  of  the  single  unit  cases.  The  out- 
put CASE  FILE  is  a tape  with  all  the  single  unit  cases  followed  by  all 
of  the  multi -unit  cases . 

MUC130  is  a cull  routine  that  extracts  all  of  the  C-130  cases 
for  either  the  east  or  west  coast  and  merges  them  at  the  end  of  the 
CASE  FILE.  This  new  file  is  called  the  PCP  CASE  FILE  and  contains 
all  of  the  single -unit  cases  followed  by  the  multi -unit  cases  and 
finally  the  C-130  cases.  Since  all  of  the  C-130  cases  are  in  the 
district  5 or  12  tape,  the  input  to  MUC130  is  either  the  district 
5 or  12  tape  and  the  CASE  FILE  from  MUTAPE. 

A sample  of  the  C-130  cases  from  the  PCP  CASE  FILE  follows: 
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III.  PCP  - The  Program  for  Calculating  Case  Parameters 

PCP,  the  computer  program  developed  to  calculate  the  case  para- 
meters required  by  the  simulation  is  comprised  of  a main  program,  PCP, 
which  calls  subroutines  NEEDS,  READ,  NUCASE  and  FIELD;  all  are  coded 
in  FORTRAN  V.  A standard  input/output  (I/O)  routine,  NTRAN1,  is  also 
called  to  allow  for  a higher  degree  of  parallel  processing,  not 
normally  permitted  by  FORTRAN  V.  By  buffering  I/O  in  parallel  with 
computing,  more  efficient  use  of  the  tape  (and/or  drum)  can  be  gained 
through  the  use  of  NTRAN. 

The  remaining  subroutines  and  main  program  descriptions  and  how 
to  use  PCP  are  discussed  in  the  succeeding  sections . Flowcharts  and  computer 
listings,  will  be  found  in  the  accompanying  NBS  Reports  10435  and  10436. 

A . The  Software 

1.  The  Main  Program  - PCP 

Most  of  the  case  parameters  are  prepared  in  the  main  program, 

PCP.  PCP  calls  subroutine  NEEDS  to  map  the  assistance  rendered 
to  the  client,  to  the  need  required  in  OPSIM  as  coded  in  the 
Resource  Capability  Matrix. 

The  inputs  to  PCP  include  the  SAR  case  descriptions  for 
the  district  undergoing  simulation,  i.e.  the  PCP  CASE  FILE 
and  the  following  data  set:  (a)  the  district  location  data;  (b) 
the  search  speeds  for  each  resource  type  in  the  inventory;  (c)  the 
needs  matrix;  (d)  the  interdistrict  (C-130)  case  location  data 
which  includes  district  centroids  and  associated  accumulative 

■HjNIVAC  1108  System  Manual , section  7,  page  9. 
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probabilities  of  occurrence  (used  for  cases  with  unknown  locations) ; 
and  (e)  the  location  of  each  OPFAC  relative  to  the  district  origin. 

A random  number  seed  is  also  required.  The  reader  is  referred  to 
Table  1,  page  28,  the  Input  Data  to  PCP^for  a description  of  each  input 
variable  delineated  above,  the  sample  input  data  which  follows  it,  and 
the  illustration  for  preparing  the  data  set  input  deck. 

Once  the  data  set  is  read  in,  PGP  processes  a case  at  a 
time  to  develop  the  simulation  required  case  parameters.  The 
reader  is  referred  to  Table  3,  page  46,  Output  Data  from  POP  for 
a description  of  each  of  these  parameters  as  they  appear  on  the 
case  tape.  Paper  output  is  also  provided  and  Figure  2,  page  48, 

Example  of  a Case  Description-Output  from  the  PCP,  illustrates 
this  form  of  output. 

Refer  to  the  Flowcharts  of  the  PCP,  Report  10435  during 
the  discussion.  A case  is  read  from  the  PCP  case  file  via  sub- 
routine NUCASE.  Some  of  the  case  data  is  converted  to  integer  for 
processing  in  PCP;  other  data  remains  in  its  alphanumeric  form. 

The  reader  should  refer  to  the  SAR  Assistance  Report  form  Figure  3, 
page  10,  during  the  discussion.  NUCASE  also  determines  CMAX,  the 
maximum  number  of  C-cards  (Reporting  Command  Data  Cases)  and  N0(I), 
the  maximum  number  of  D-cards  (Assisting-Resource  Data)  associated 
with  each  C-card. 

The  multi-unit  case  control  number,  A3,  is  converted  to 
integer;  the  case  number  is  retained  and  therefore  only  the 
alphanumeric  is  changed.  This  conversion  is  listed  below: 
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CONVERSION  OF  A3  TO  INTEGER 


Old  A3 

New  A3 

DXXXX 

1XXXX 

UXXXX 

2XXXX 

JXXXX 

3XXXX 

KXXXX 

4XXXX 

LXXXX 

5XXXX 

MX  XXX 

6XXXX 

NXXXX 

7XXXX 

OXXXX 

8XXXX 

PXXXX 

9XXXX 

QXXXX 

11 XXX 

RXXXX 

12XXX 

SXXXX 

13XXX 

TXXXX 

I4XXX 

EXXXX 

15XXX 

WXXXX 

16XXX 
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Unit  and  C-130  Case  Format 


. I Ki  AbDK’Y  I’l  I' Al  i I • I-  N T 
I u-  s.  COAST  GUARD 
CG-327  2 (Rev.  9-65) 


! REPORT 


ASSISTANCE  REPORT 

M s s,  v c\ r \a »\v  or-  R-rc  l c gy ct u r u. ld 


• REPORTING  UNIT  (In  cl  ude  lo  c al  ion  o r p amion  on  t ol  ol  i onj 

j _ <5030  ■c>>' 7 „ Cx  Co  u c &?<**,  /VjgAd- 


Lf_  Col,  ] a. 

1 |0 


IDENTIFYING  DAT Ai 


Xb 


11 


1 I Zj_L06ICiH_  fields  plc,  CnOE 

o.  / 


O „ / ,3  ,L 


Unit  Accounting  Code 


O 


U 


^ ^ j Unit  Case 


Number 


o 2.  , 2 V 


(9, 


o 


Multi-Unit  Case  Controller  Code 
and  Case  Number 


Month  and  Yenr  Unit  Notified  (Item  Cl) 


instruction:; 

Reference:  Latest  revision  of  COMDTINST  312 3.9 

1.  Prepare  original  only.  Use  of  pen  or  pencil  in  permitted. 

2.  Shaded  areas  require  uso  of  code  list  at  oil  times. 

3b  ALL  SPACES  MUSI'  BE  I-ILLED.  See  reference  for  codes 
when  an  item  is  unknown  or  not  applicable. 

* Y - yes,  N - no,  U - Unknown  or  not  applicable;  DIR  - Defi- 
ciency Improvement  Report 
**  To  the  nearest  tenth  of  on  hour. 


0.  DISTRESSED  UNIT  DATA:  CARD  NO.  1 


21 

1 [21 

0 , /I 

s//.3.a,;  . 

A 

Distressed  Craft 
Number 

XL  12 

62 

o\ 

it 

P 

Type  of  Di stress 

21 

2 

29 

I \ o 

— - 1 Date/Time  Distress  Incident 
& , [Occurred 

<■*6  13 

65 

o 

u 

Dis 

tressed  Unit  Description 

25 

3 

35 

op 

Total  Number  of  Persons  in  Distress  ' 
or  on  Distressed  Croft 

C'? 

14 

67 

S' 

Source  of  Information 

Total  Number  of  Lives  Lost  Among 

62 

V\ 

4 

39 

L ^ i l 1 

Persons  in  Distress 

15 

68 

o 

Means  of  Notifying  Const  Guard 

5 

43 

o,  o 

OJ_ 

S.o  ,o 

0 

Value  of  Property 
Involved 

70 

16 

70 

o 

o 

7 

Type  of  Distress  Area 

51 

6 

51 

S,  o 

Air  Temperature  (°F) 

Record  the  severest 

T6 

17 

73 

o 

Nature  and  Severity 

NO 

_n 

7 

53 

'fp 

Sea  Temperature  (°F) 

weather  encountered  in  the 
distressed  unit's  area. 
(Required  in  all  coses 
unless  omission  is  specif- 
ically covered  in  the 
reference.) 

NA  - not  applicable 

18 

76 

1 

Cause  ol  Incident 

i r 

8 

55 

2,o 

Wind  Force  (ICpots) 

17 

19 

77 

3 

Final  Status  of  Property 

: ST 

9 

57 

o y Co 

Height  of  Sea  (Feet) 

7'i 

20 

78 

r 

Was  Distressed  Unit’s  Emergency  Equipment 
Satisfactory?  (No  - COMMENT'  REQUIRED) 

si 

10 

52. 

0.3 

Visibil  ity 
_(Ncarcst_Milc)  . 

UK 

- unknown 

7^ 

21 

79 

>\f 

Were  Non-Coast  Guard  Units  in  a Position  and 
Capable  of  Performing  Rescue? 

i fcl 

l 1 

61 

\ j 1 Was  Distressed  Unit  Being  Set  Onto  a Dangerous 
/ j Bench  or  Shoal  Within  5 Miles? 

10 

22 

10 

PL 

Vessels  65  Feet  and  Under;  Bad  Vessel 
been  Boarded  within  Past  Year? 

C.  REPORTING  COMMAND  DATA:  Card  No.  2 (Cover  ONLY  Reporting  Command  Porformonco  except  Items  C6  on<J  C7.) 

ThiVs  ui stressc ti  Unit” io w e <I7" Ls co rte  d 


1 21 


/ v5"  / ,<9,  / ,5“ 


Date/Time  Notified 


lit 


i mui  ) i * 

9 ,56  O i O t CP  O j or  Transported 


'll 


11 

_ CV3 

102 


CJ  Day  of  Week  When  Notified, 


iao 


O, 


Number  of  Bodies  Recovered 


Number  of  Assistance  Cases  in  Progress  at  Time 
Notified 


Ui 


11  63 


Number  of  Persons  Saved 


/ 


'Date/Time  First  Assisting 
Resource  Underway 


12|  66 


O 


Number  of  Persons  Otherwise  Assisted 


Number  of  Assisting  Craft  or  Groups  Assigned  to 
Case 


lZf? 


13  69 


Number  of  Persons  Indirectly  Assisted 


36 


- Dutc/'i  ime  l'  irst  Assistance 
/ , O / , CJ  , -c  C>  at  Distress  Position 


l7- 


/ ,-5~  1 / , / ,3s 


Dnte/Time  Case  Terminated 


10 


113 


12 


. LATITUDE  . 

V,  / I3.VIV 


LUNG 

<9,7  ,Q 


LONGITUDE 

7 ,0 


0,0  o 


Position  nt  which  Distressed  , 

Unit  Located.  (Report  even  it  j 

Locati  d by  another  Unit.)  j | a M 


78 


0 I Report  Review  Requirements 


Distressed  Unit’s  Distance  (miles)  from 
Reported  Position 


14  0 


iU 

113 


1 11 T 

lb  j 
1,61 
160 
lot 
rn 

IT  4 


D.  ASSISTING  RESOURCE  DATA:  Cord  No.  3 (Ono  card  for  each  rosourco.) 

1 2lL  I ,3  j ’Vype  of  Resource  Reported  in  this  Section 


ISO 


23  O , V 


33  / ,-5~ 


V .3  ,7  ,7 


/ o / p 


/.O  2.  ,5” 


0 0 0 2 


. 8 54 


o 


Assisting  Boat,  Vessel  or 
Aircraft  Number  


Dnto/Timo  (3)  Undorwny, 

(4)  On  Scene,  Slnrtcd  Searching 
or  Rendezvoused. 


Distance  to  Scene,  Search  Aroa  or 
Rendezvous  Point 


142 


It  4 


1'iG 

m 


Thi.  item  for  HQ  USE  ONLY 


02 


14 


Hours  Spent  Searching  or  Conducting  , o 

Communications  Check  


O Q OO  (a 


Total  Miles  Traveled  on  Cose 


91  58  O ■ ! 


Q O / . 3-  Total  Time  Underway  on  Case 


Number  of  Sorties 


15 


IS  2. 


OJ 


U- 


68 


70 


72 


73 


O 2 


/ 0 


L 


Method  of  Locating,  Distressed  Unit 


Assistance  Rendered  to  Personnel 
Prlumiy  Annin limcn  Rendered  Relniing  to 

Property . 

"Secondary  Assistance  Rendered  Relating 
to  Prop  City 


On  Scene  Frcqviency 


Ship/ Shore  or  Air/ Surface  Control  Circuit 


Was  Equipment  and/or  Personnel  Performance 
Satisfactory?  (No  - DIR  REQUIRED) 
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Figure  3 


The  location  of  the  cases  (x,y)  relative  to  the  district 
origin  (LAT,  LON)  is  determined  next.  If  the  location  data 
C7A(I)  or  C7B(I)  is  a zero  or  9999  and  no  location  data  is  given 
for  the  OPFAC  A1B(X) ; i.e.,  it  is  undefined  in  SARSIM,  then  the  case  is 
assigned  the  following  coordinates  (XPT,  YPT)  relative  to  the 

district  origin  (note:  these  can  never  be  (0,0):  The  above  is 

true  for  both  single  unit  and  multi-unit  cases  (in  single  unit  cases,  1=1.) 

If,  however,  at  least  one  C-card  does  contain  the  location  data, 

that  location  is  used  relative  to  the  district  origin,  (LAT,  LONG) . 

If  the  case  does  not  have  location  data  but  does  have  a defined 
OPFAC  A1B(I)  with  a location,  the  case  then  takes  on  the  location 
(XDELTA,  YDELTA) , from  A1B(I) . 

If  the  case  does  have  a valid  location  (C7A(I) , C7B(I))  then 
(x,y)  is  calculated  using  plane  geometry. 

If  the  case  is  an  interdistrict  (C-130)  case,  with  a valid 
(C7A(I) , C7B(I))  then  the  distance  R to  each  district  centroid, 
is  taken  as  the  minimum  R.  Once  A1A(I)  has  been  determined,  (x,y) 
is  calculated.  If  the  case  occurred  in  the  district  being 
exercised,  DIST,  then  A1B(I)  is  set  to  zero.  (A  zero  code  in 
A1B(I)  causes  OPSIM  to  assign  the  case  to  the  closer  station  in 
DIST.) 

If  an  invalid  (C7A(I) , C7B(I))  is  given  for  interdistrict 
C-130  cases,  then  A1A(I)  is  assigned  via  the  Monte  Carlo  technique 
from  the  accumulative  probabilities,  PROB(D) , where  D is  the 
district . 
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Once  (x,y)  has  been  found  for  all  cases  in  DIST,  a check 
is  made  to  see  of  (x,y)  is  in  the  range  (XLOW,  XLHT)  and  (YLOW, 

YLMT  . If  not,  (x,y)  is  arbitrarily  set  to  (2,2). 

Note,  that  for  multi -unit  cases  with  multiple  C- cards  the 
(x,y)  is  retained  in  (xx,yy)  until  all  C-cards  are  reviewed. 

MIMCI  is  the  date  and  time  of  notification  for  the  case  and 
is  taken  as  Cl(l),  the  date  and  time  reporting  command  was  notified, 
of  the  first  C-card.  (The  cases  are  sorted  chronologically  within 
each  set  of  SAR  assistance  forms  for  each  multi-unit  case.) 

STANO,  the  sequential  number  of  the  primary  station,  becomes 
A1B(1)  the  OPFAC  corresponding  to  Cl(l).  If  A1B(1)  is  not  in 
the  set  of  current  existing  OPFAC' s,  then  a new  STANO  is  created, 
added  to  the  list  and  STANO  becomes  the  previous  STANO  number  plus 
one.  If  no  A1B(1)  exists,  i.e.  it  is  zero,  then  STANO  is  zero. 

Next  TSM,  the  total  number  of  search  miles  expended  on  the 
case  is  found  as,  TSM  = D"c*rds  (D6(I,J)  *S0A3  (D1(I,J))  where 
D6(I,J)  is  the  hours  spent  searching  (by  the  Ith  resource  from 
the  Jth  Opfac)  and  S0A3  (D1(I,J))  is  the  search  speed  of  the 
assisting  resource  type  Dl(I,J). 

A simple  conversion  is  made  of  B16,  the  code  for  location 
of  the  case  in  terms  of  miles  off-shore  is  made.  Called  MILES, 
the  conversion  is  given  below: 


B16  Code 

MILES  (tenths  of  nautical  miles) 

B16A  = SXX 

0.0 

> 1,  < 6 

3.0 

= 7,  = 8 

30.0 
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80.0 


= 9 

= 10,  . ..,  998  actual  number  of  miles 
* 999  9.0 

Next,  B17B,  the  coded  severity  of  the  case  is  translated 


thus,  to  NSEV: 

B17B  Code  NSEV 

1 1 

2 2 

3 4 

4 3 

5 5 


If  the  total  search  time  D6T0T,  calculated  over  all  the 
D-cards  on  the  case  exceeds  five  tenths  of  an  hour,  the  case  is 
considered  a short  search  case;  SI  is  set  to  zero  and  if  the 
first  resource  searched  (first  D6(I,J)  is  greater  than  zero)  then 
S2  is  set  to  one.  If  any  other  resource  other  than  the  first 
searched,  then  S2  is  set  to  two. 

For  each  resource,  the  time  the  resource  arrives  on  scene 
and  commences  serving  a need,  B4FRST  after  searching  is  found: 
D4FRST(I ,J)  * D4(I,J)  + D6(I,J),  where  D4(I,J)  is  the 
date  and  time  the  resource  arrived  on  scene,  and  D6(I,J)  is  the 
time  spent  searching.  Next  subroutine  NEEDS  is  called  to  trans- 
late the  assistance  rendered  to  the  proper  NEED(N)  code.  If  the 
need  of  a case  NEED(N)  is  any  of  the  following:  10,  general 
escort;  15,  tow;  17,  air  escort;  or  19,  rescue  and  tow,  then  the 
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number  of  needs  other  than  tow  N is  decremented,  the  number  of 
tow  resources  M is  incremented  and  the  need  stored  in  MEED(M) . 

(Any  miscellaneous -classed  assisting  resource  data,  i.e.,  any 
D1(I,J)  > 70,  is  not  included  in  the  output  case  description). 

The  array  NEED(N)  is  a set  of  requirements  for  aid  other 
than  tow  for  each  case.  The  array  MEED(M)  is  a temporary  array 
of  tow  (including  escort)  needs.  If  a case  has  a requirement  for 
tow  (M>0) , only  one  such  need  can  exist.  Therefore,  the  MEED(M) 
array  is  searched  for  the  need  with  the  highest  rank  and  NEED(NN) , 

NN  = N+l,  is  set  to  this  need.  The  ranking  order  from  high  to 
low  is:  17,  15,  19,  10. 

T0S(N) , the  time  spent  on  scene  serving  NEED(N)  is  calculated: 
TOS(N)  = D8(I , J)  - D6(I,J)  - ENRTE  *2*D9(I,J);  where  ENRTE  = 

D4(I,J)  - D3(I,J)  and  D8(I,J)  = Total  time  underway  for  this  resource. 

D6(I,J)  = Hours  spent  searching  for  this  resource. 

D9(I,J)  = Number  of  sorties  for  this  resource. 

D4(I,J)  = Date/Time  on  scene  to  start  searching  for  this 

resource . 

D3(I,J)  = Date/Time  underway  for  this  resource. 

Finally,  for  this  assisting  resource,  the  time  the  resource  leaves 
the  scene  after  performing  his  need,  LVTIME  (I,J)  is  calculated. 

The  sets  D4FRST(I,J)  and  LVTIME (I ,J) are  used  to  find  TE, 

MAX 

the  elapsed  time  of  the  case.  LAST  is  the  j j (LVTIME  (I,J)} 
and  FRST  is  the  (D4FRST  (I,J)}.  The  elapsed  time,  TE,  is 

given  by: 

TE  = LAST  - FRST 
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TE  is  used  to  calculate  DELTA(N) , for  each  need  of  the  case , 

The  MAX  (T0S(N)},  MAXTQS,  and  the  E TOS(N),  SUMTOS (N) 

N 

are  calculated,  NN  is  set  and  GAM4A,  the  degree  of  non-parallel 

» -r 

service  on  a multi -resource  case  is  calculated,  Prior  to  output  of 
the  case,  checks  are  made  on  the  data  for  such  crucial  parameters 
as  M,  N,  and  associated  NEED(N)  and  MEED(M)  • In  short,  general 
housekeeping  is  performed  on  each  case  such  that  no  illegal  in- 
formation is  left  to  cause  problems  when  simulating  the  case  in 
OPSIM. 

So  far  the  discussion  has  stemmed  from  processing  cases  with 
no  search  or  short  search  involved.  For  long  search  cases,  (i.e., 
D6T0T,  the  total  hours  in  search,  greater  than  or  equal  to  .5) 
special  processing  is  required  to  calculate  the  parameters 
describing  long  searches.  If  the  case  is  a single -resource  case, 
(i.e.,  CMAX  and  NO(X)  are  both  one)  and  D6T0T  does  not  exceed  14 
hours,  then  SI,  the  number  of  search  resources  is  set  to  one.  If 
on  the  other  hand,  this  is  a multi -resource  long  search,  (i.e.’, 
either  CMAX  qr  NO (I)  exceed  one),  and  D6T0T  does  not  exceed  14 
hours,  NITEHR,  the  nighthours  elapsed  during  the  case,  is  set  to 
zero  and  SI  calculated:  SI  = D6T0T/ (elapsed  time  - NITEHR).  If 

D6T0T  exceeds  14  hours,  NITEHR,  is  calculated  and  subtracted  from 
elapsed  time  to  determine  SI,  (A  constant  0.99  is  added  to  SI  so 
that  SI  becomes  the  next  integer  value  after  truncation) . For 
all  long  search  cases,  a flag  S2  is  set  to  0. 


-15- 


One  additional  conversion  is  performed  on  the  parameter 
B13,  the  client  description,  which  contains  the  coded  length 
of  the  vessel.  These  conversions  are  listed  below: 


Old  B13 

New  B13 

50 

66 

51 

101 

52 

201 

60-99 

0 

The  remaining  parameters  B3,  B5,  B6,  B8,  B9,  BIO,  and  B12A 
are  ascertained  in  Subroutine  FIELD,  and  not  changed  in  PCP. 
See  page  24  for  an  explanation  of  these  parameters. 

Once  the  EOF  is  reached,  the  PCP  halts. 
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2)  Subroutine  NEEDS 


Sib rout in©  MEEDS  performs  a translation  of  the  historical 
assistance  r«tel®red  (as  reported  in  the  SAR  assistance  form)  to  the  capa- 
bilities delineated  in  the  needs  portion  of  the  Resource  Capability  Matrix. 

(The  reader  is  referred  to  the  analyst  documentation  for  a review  of 
this  matrix. ) 

Generally,  the  assistance  rendered  to  personnel,  for  each  resource 
serving  a case  is  queried  before  the  primary  and  secondary  assistance 
rendered  to  property,  for  the  translation.  This  query ^coup led  with  an 
investigation  of  the  historical  resource  rendering  aid,  forms  the  core 
of  NEEDS.  INEED(N)  is  the  input  required  by  NEEDS.  The  reader  is  referred 
to  Table  1 of  the  Input  Data  to  PCP,  the  NEEDS  matrix,  for  the  correlation  of 
Dll  and  D12  between  the  code  and  the  Resource  Capability  Matrix.  Each 
time  a D-card  (assisting  resource  data  card)  for  a case  is  encountered, 
in  the  main  program,  PCP,  subroutine  NEEDS  is  called.  Upon  entry  to  NEEDS, 
the  assisting  resource  code  is  queried  to  see  if  an  aircraft  rendered 
assistance^  If  so,  and  if  the  aircraft  assisted  personnel  (i.e.,  D11(I,J)  > 0) 
then  a need  of  the  case  being  processed,  NEED(N)  becomes  INEED(D)  , where  D 
= D11(T,J),  If  the  aircraft  instead  rendered  primary  assistance  to 
property  (D12(I,J)  > 0),  then  NEED(N)  becomes  INEED(D)  , where  D = D12(I,J). 

If  the  client  description,  B13,  indicates  an  aircraft  and  D12(I,J)  is 
coded  68,  escort,  then  the  NEED(N)  is  set  to  17,  indicating  the  client  re- 
quires an  escort  by  an  aircraft.  See  Table  2,  page  18.  Explanation  of 
NEED(N)  for  code  description.  If  an  aircraft  was  called  to  scene,  but  no 
assistance  rendered,  i.e.,  Dll  and  D12  are  zero,  then  the  NEED  is  set 
to  18. 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 


Explanation  of  NEED(N) 


Description 
Provide  Equipment 
Deliver  Pump/Equipment 
Made  Repairs 
Fought  Fire 
Vector  Other  Unit 
Dewater 
Refloat 
Icebreaking 
Refueled  and  Supplied 
Surface  Escort 
Standby 

Located  Property  and  Owners  Advised 
Freed  from  Position  of  Peril 
General  Assistance  Rendered  Surface 
Tow;  Dev/at er  and  Tow;  Refloated  and  To’ 
Evacuate  (Rescue)  People  on  Board 
Air  Escort 

General  Assistance  Rendered 
Rescue  and  Towed 


Table  2 
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If  the  assisting  resource  was  a surface  resource,  Dll  > 0,  and  either 
D12  or  D13  were  one  of  the  following;  tow,  dewater  and  tow,  refloated 
and  tow,  tugbird  tow,  or  recovered  property  and  returned  to  owner,  then 
NEED(N)  is  set  to  19  implying  that  assistance  was  rendered  to  personnel 
and  property  as  described  above. 

If  the  surface  resource  assisted  property  only,  and  D12  was  any  one 
of  the  following: tow,  dewater  and  tow,  refloated  and  tow,  tugbird  tow 
or  recovered  property  and  returned  to  owner,  then  NEED(N)  is  set  to  15, 
a requirement  for  tow.  If  however  D12  was  none  of  these  described,  then 
the  NEED(N)  is  extracted  from  the  INEED(D)  matrix  unless  the  assistance 
rendered  was  escort  and  the  client  an  aircraft.  In  the  latter  case,  NEED(N) 
is  set  to  17. 

If  the  surface  resource  is  called  to  scene,  and  no  assistance  is 
rendered  to  either  personnel  or  property  then  NEED(N)  is  set  to  14,  or 
general  assistance  rendered  by  a surface  resource. 

Once  NEED(N)  is  set,  control  is  passed  back  to  the  calling  main 
program,  PCP. 
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3)  Subroutine  NUCASE 


Unit  and  C130  cases  are  processed  in  groups  of  32  words  (logical 
fields,  LF's).  The  first  logical  field  is  comprised  of  the  number  of 
LF's  in  the  case,  the  A-card,  B-card,  C-card,  and  one  D-card.  The  next 
logical  field  contains  the  number  of  D-cards  and  up  to  3 D-cards.  Each 
successive  LF  contains  up  to  3 additional  D-cards.  If  a case  reports 
no  D-cards  the  number  of  LF's  in  the  case  is  set  to  zero. 

Multi-unit  (MU)  cases  have  an  extra  logical  field  at  the  beginning 
of  the  case  which  gives  the  number  of  LF’s  in  the  case,  the  number  of 
participants  (C-cards)  and  the  B-card.  The  number  of  LF's  for  each 
participant  with  the  corresponding  A-card,  C-card  and  D-card  data  follows, 
the  participants  sorted  by  date  and  time  of  notification. 

This  routine  calls  READ  to  get  a buffer  of  data  (LOAD)  ready  for 
processing.  LP  is  then  used  as  a pointer  to  indicate  the  word  in  LOAD 
where  the  case  begins  (LOAD(LP)).  Since  a case  may  exist  between  buffers, 
it  is  transferred  to  an  array  called  DATA  for  further  processing. 

NUCASE  then  calls  FIELD  to  extract  the  data  to  be  used  by  the  main 
program,  PCP.  The  number  of  C-cards  (CMAX)  and  the  number  of  D-cards  for 
each  C-card  (NO (I))  is  set  for  each  case. 

Figure  4 relates  the  first  logical  fields  and  successive  fields  to 
the  SAR  assistance  form. 
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TREASURY;  DEPARTMENT 
j.  & COAST  GUARD 
1G-327  2 (Rev.  9-65) 


MS^ 


ASSISTANCE  REPORT 
FORMAT  OP  FIRST  LOSICRU  FUlt 


j 

: REPORT  CONTROL  SYMnoi. 

osfi^ooc 

i 

ALU  “TIMES  LOCAL 


REPORTING  UNIT  (Include  location  or  permanent  station) 


Cr  /&  gc^st; e,iT  S£e5  O o O , & J o u c <as  fr/sr , 


F COL. 

W — — - 1 ^ . - J 

A.  IDENTIFYING  DATA:  | 

INSTRUCTIONS 

1 

0 

1 

2 

LOSICftL  FIllSS  PER.  ChSE 

Reference:  Lntcst  revision  of  COMDTINST  3123.9 
L.  Prepare  original  only.  Use  of  pen  or  pencil  is  permitted. 

2.  Shaded  areas  require  use  of  code  list  at  all  times. 

3.  ALL  SPACES  MUST  BE  FILLED.  See  reference  for  codes  ; 
when  an  item  is  unknown  or  not  applicable. 

* Y • yes,  N » no,  U - Unknown  or  not  applicable;  DIR  - Defi- 
ciency Improvement  Report 
**  To  the  nearest  tenth  of  an  hour. 

2 

1 

2 

OJ_. 

3,0 

/ .3  .6  1 

Unit  Accounting  Code 

9 

2 

9 

O 

2.2 -A. 

Unit  Ca’fee  Number 

15 

3 

13 

u 

o 

2.  , 2 , V 

Multi-Unit  Case  Controller  Code 
and  Case  Number 

n 

4 

18 

o, 

3 

(o 

Month  ai 

id  Year  Unit  Notified  (Item  Cl) 

B.  DISTRESSED-UNIT  DATA:  CARD  NO.  1 

21 

1 

21 

oxn 

S ^,3,2 

/ .A 

Distressed  Craft  i 

Number 

42  j 

12 

52 

0.4 

/ 

1 

Type  of  Distress 

21 

2 

29 

1 ,5 

4>,5“ 

Date/ Time  Distress  Incident 
Occurred 

45> 

13 

55 

2 

H 

Distressed  Unit  Description 

IS 

3 

35 

O',o  ,o  ,S 

Total  Number  of  Persons  in  Distress  ' 
or  on  Distressed  Craft 

41 

14 

57 

S 

Source  of  Information 

11 

4 

39 

O,—, — . 

Total  Number  of  Lives  Lost  Among 
Persons  in  Distress 

(o% 

15 

68 

0,-2 L 

Means  of  Notifying  Coast  Guard 

AS 

5 

43 

0,0,0,  / ,£~ ,0,0 

o 

Value  of  Property 
Involved 

1 

70  16 

70 

o.o  .9 

i 

Type  of  Distress  Area  ; 

51 

6 

51 

Sp 

Air  Temperature  (°F) 

Record  the  severest 

15 

17 

73 

o 

6 ,s 

Nature  and  Severity 

55 

7 

53 

4p 

Sea  Temperature  (°F) 

weather  encountered  in  the 
distressed  unit’s  area. 
(Required  in  all  cases 
unless  omission  is  specif- 
ically covered  in  the 
reference.) 

NA  - not  applicable 
UK  - unknown 

lb 

18 

76 

1 

Cause  of  Incident 

55 

8 

55 

2,0 

Wind  Force  (Knots) 

11 

19 

77 

3 

. f 

Final  Status  of  Property  | 

51 

9 

57 

O , <£> 

Height  of  Sea  (Feet) 

-n 

20 

78 

r 

Was  Distressed  Unit’s  Emergency  Equipment 
Satisfactory?  (No  - COMMENT  REQUIRED) 



fcl 

:o 

59 

3 

Visibility 
1 N earest  Mile) 

19 

21 

79 

a r 

Were  Non-Coast  Guard  Units  in  a Position  and 
Capable  of  Performing  Rescue? 

li 

61 

* 

Was  Distressed  Unit  Being  Set  Onto  a Dangerous 
Beach  or  Shoal  Within  5 Miles? 

SO  j 22 

80 

y 

Vessels  65  Feet  and  Under:  Had  Vessel 

been  Boarded  within  Past  Year0 

C.  REPORTING  COMMAND  DATA:  Card  No.  2 (Cover  ONLY  Reporting  Command  Performance  except  Items  C6  and  C7„) 

21 

i 

21 

/ 

3T 

1,0,/ ,s 

Date/Time  Notified  ' 12(> 

9 

56 

0,0,0 

Miles  Distressed  Unit  Towed,  Escorted  | 
3 or  Transported 

n 

2 

27 

3 

Day  of  W’eck  When  Notified 

10 

60 

O,  — 

Number  of  Bodies  Recovered 

n 

3 

28 

6? 

Number  of  Assistance  Cases  in  Progress  at  Time 
Notified  ,12  5 

11 

63 

0 ,o  ,S 

Number  of  Persons  Saved 

29 

4 

29 

/ 

& 

/ ,o,/ ,s 

■Date/Time  First  Assisting 
Resource  Underway 

124 

12 

66 

or— 

Number  of  Persons  Otherwise  Assisted 

95 

5 

35 

/ 

Number  of  Assisting  Craft  or  Groups  Assigned  to 
Case 

121 

13 

69 

0,-7- 

Number  of  Persons  Indirectly  Assisted 

it 

6 

36 

/ .s 

/ ,0,2. 

jr 

Date/Time  First  Assistance 
at  Distress  Position 

152 

14 

72 

i 

/ , / ,3  ,«S  Date/Time  Case  Terminated 

102 

7 

4? 

, . L ATITUDE  , ^ f 

4,  / 3,V  V 

Position  at  which  Distressed 
Unit  Located.  (Report  even  if 
Located  by  another  Unit.) 

m 

15 

78 

o 

Report  Review  Requirements 

101 

47 

LONGI 

0,7  .O 

TUDE 

4 ,o\hj 

151 

16 

79 

This.  Item  for  HQ  USE  ONLY 

115 

8 

53 

0.0,0 

Distressed  Unit’s  Distance  (miles)  from 
Reported  Position 

140 

117 

80 

D. 

A 

SSISTING  R ESC 

)URCE  DATA:  Cord  No.  3 (On©  card  for  each  rosourcc.) 

• 

1 

1 

?1 

7 3 

1 Tvoe  of  Resource  Reported  in  this  Section 

140  To 

60 

O,  ! 

Method  of  Locating  Distressed  Unit 

M3 

2 

O' 

V 

- V ,3  , 7 ,7 

Assisting  Boat,  Vessel  or 
Aircraft  Number 

1«2  11 

62 

/,/ 

Assistance  Rendered  to  Personnel 

111 

3 

2S 

ZZ-Jl i— 

/,<s- 

* ,0J  f 

Date/Time  (3)  Underway, 

134  4.2 

64 

(o  p 

Primary  Assistance  Rendered  Relating  to 
Property 

(4) 

On  Scene,  Started  searemng 
Rendezvoused. 

14b 

13 

66 

op 

Secondary  Assistance  Rendered  Relating 
to  Property 

155 

4 

V 

/ 1*5* 

/ o o.S 

or 

— 

Ibl 

5 

4 

£ 

«0  o 

o z 

Distanc 

Rendez 

e to  Scene,  Search  Area  or 
vous  Point 

in 

14 

68 

O 2 
> — 

On  Scene  Frequency 

IfcS 

1 L- 1 k. 

Hours  Spent  Searching  or  Conducting 

Ho 

15 

70 

/ 

o 

Ship/Shore  or  Air/Surface  Control  Circuit 

1W 

7 

4 

1 

5 O <D  (?,0 

Total  Miles  Traveled  on  Case 

112 

16 

72 

t 

Was  Equipment  and/or  Personnel  reriormante 
Satisfactory?  (No  • DIR  REQUIRED) 

— 

US 

8 

. — ■ — k_ 

On  1 7. 

To 

tal  Time  Underway  on  Case 

17 

73 

114 

9 

5 

8 0,  / 

Number 

of  Sorties  _ ...  - „ 

18 

\ 

•-  ' 

Fimas'  4* 
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NB5 

format  for  successive  logical  ptclks 

1 1 

0 

NUM&Efc  Of  3CARD5  PER.  CftSt 

3 

1 

21 

Type  of  Resource  Reported  in  this  Section 

42  ! 

10 

30 

Method  of  Locating  Distressed  Unit 

i 5 

2 

23 

i 

Assisting  Boat,  Vessel  or 
Aircraft  Number 

44 

11 

62 

1 

Assistance  Rendered  to  Personnel 

I 11 

3 

29 

i 

Date/Time  (3)  Underway, 

(4)  On  Scene,  Started  Searching 
or  Rendezvoused 

4fc 

12 

64 

<, 

Primary  Assistance  Rendered  Relating  to 
Property 

17 

4 

35 

4% 

13 

66 

1 

Secondary  Assistance  Rendered  Relating 
to  Property 

Z*> 

5 

41 

Distance 

Rendezv 

to  Scene,  Search  Area  or 
ous  Point 

50 

14 

68 

1 

On  Scene  Frequency 

21 

6 

45 

** 

Hours  Spent  Searching  or  Conducting 
Communications  Check 

52 

15 

70 

Ship/Shore  or  Air/ Surface  Control  Circuit 

31 

7 

49 

‘ 

Total  Miles  Traveled  on  Case 

54 

16 

72 

* 

Was  Equipment  and/or  Personnel  Performance 
Satisfactory?  (No  • DIR  REQUIRED) 

6fa 

8 

54 

L~ — 1 “-T»- 

Total  Time  Underway  on  Case 

55 

17 

73 

40 

9 

58 

Number  of  Sorties 

5fc> 

18 

D. 

ASSISTING  RESOURCE  DATA:  Cord  No.  3 (On.  card  lor  each  roiourco.)  | D.  ASSISTING  RESOURCE  OAT Al  Cord  No.  3 (On.  cord  lor  ooch  r.lourc.l 

i 

! 

1 

21 

* 

Type  Of  Resource  Reported  In  this  Section 

10 

60 

1 

Method  of  Locating  Distressed  Unit 

1 

i 

*2 

23 

Assisting  Boat,  Vessel  or 
Aircraft  Number 

4% 

11 

62 

Assistance  Rendered  to  Personnel 

(o5 

3 

29 

Date/Time  (3)  Underway, 

(4)  On  Scene,  Started  Searching 
or  Rendezvoused 

100 

12 

64 

Primary  Assistance  Rendered  Relating  to 
Property  j 

It 

4 

35 

10  2 

13 

66 

Secondary  Assistance  Rendered  Relating 
to  Property 

77 

5 

41 

I I 1 " 

Distance  to  Scene,  Search  Area  or 
Rendezvous  Point 

10H 

14 

68 

On  Scene  Frequency 

6 

45 

e* 

Hours  Spent  Searching  or  Conducting 
Communications  Check 

lCIo 

15 

70 

Ship/Shore  or  Air/ Surface  Control  Circuit 

! vs 

7 

49 

Total  Miles  Traveled  on  Case 

los 

16 

72 

* 

Was  Equipment  and/or  Personnel  Performance 
Satisfactory?  (No  - DIR  REQUIRED) 

"10 

’t 

34 

** 

Total  Time  Underway  on  Case 

104 

17 

73 

1 -M 

.58 

Number  of  Sorties 

110 

18 

D.  ASSISTING  RESOURCE  DATA:  Cord  No.  3 (One  cord  for  each  resource.) 

Ui 

i 

21 

i 

Type  of  Resource  Reported  in  this  Section 

150 

10 

60 

1 

Method  of  Locating  Distressed  Unit 

113 

2 

23 

J 

Assisting  Boat,  Vessel  or 
Aircraft  Number 

15.2 

11 

62 

1 

Assistance  Rendered  to  Personnel 

m 

3 

29 

Date/Time  (3)  Underway, 

154 

12 

64 

Primary  Assistance  Rendered  Relating  to 
Property 

125 

4 

35 

i 

or  Rendezvoused 

15fc 

13 

66 

Secondary  Assistance  Rendered  Relating 
to  Property 

131 

5 

41 

- i-  i . l * 

Distance  to  Scene,  Search  Area  or 
Rendezvous  Point 

15S 

14 

68 

On  Scene  Frequency 

13.5 

6 

45 

#« 

Hours  Spent  Searching  or  Conducting 
Communications  Check 

IfeO 

15 

70 

Ship/ Shore  or  Air/ Surface  Control  Circuit 

154 

7 

49 

Total  Milea  Traveled  on  Case 

K2 

16 

72 

* 

Was  Equipment  and/or  Personnel  Performance 
Satisfactory?  (No  - DIR  REQUIRED) 

144 

y 

53 

*« 

Total  Time  Underway  on  Case 

H*> 

17 

72 

• 

14* 

y 

58 

! 

Number  of  Sorties 

H4 

18 

HAVE  YOU  PREPARED  ONE  SECTION  D FOR  EACH  CRAFT  OR  GROUP  INDICATED  IN  ITEM  CS? 

NAME  OF  DISTRESSED  CRAFT 

Alsry  /r^y 

NAME  AND  ADDRESS  OF  OPERATOR 

A.  0.  See 

/O’O.o  A? a •**  Str<e^'C 
M 9r&  . 

NAME  AND  ADDRESS  OF  OWNE«/lF  DIFFERENT 

COMMSNTO  (Llot  namota  md  addraizaao  ot  persona  reported  In  CIO  and  CJJ ) 

&•  ^6 S “ Same  &s  a 

J'os&.fh  Seg  — * ~ 

$•  S&B/ Jr.  “ 

S&GH  — 

j:  jiwss  - /^aa.’  <5*‘,  &WflSer 


: 

DATE  SUBMITTED 

<810  NATURE 

OAT* 

INITIALS 

| 

3/tC/g6 

A6&"  l.  S£~4Ma/ 

REVIEWED  AND 
FORWARDED 
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4.  Subroutine  READ 

Records  are  read  from  the  input  tape,  PCP  Case  File,  in  NTRAN 
which  allows  for  buffered  I/O  in  parallel  with  computation.  256 
words  are  read  into  one  buffer  while  the  other  buffer  is  being  processed. 
Refer  to  the  READ  listing  in  Report  10436. 

NPREC  indicates  the  number  of  records  (256  words)  that  have  been 
read  for  a particular  OPFAC.  If  it  is  set  to  zero  in  NUCASE  the  relevant 
information  from  from  a new  OPFAC  header  (32  words)  must  be  saved,  i.e. 
IPREC  set  to  the  number  of  physical  records  in  that  OPFAC.  NPREC  is 
incremented  every  time  a new  record  is  read,  and  when  it  becomes  equal 
to  IPREC  the  OPFAC  is  completed  and  NPREC  set  backs  to  zero  in  NUCASE. 

When  all  the  OPFAC  data  has  been  read  and  processed  a MU  header 
will  be  encountered.  In  addition  to  setting  IPREC,  a MU  flag  is  set 
to  let  the  main  program  know  that  MU  cases  are  being  processed  (MFLAG=1) . 
After  MU  cases  have  been  completed  the  C-130  header  is  encountered  and 
a C-130  flag  set  (CFLAG=1) . 

When  an  EOF  is  reached,  this  routine  returns  to  a specified  state- 
ment number  in  the  main  program. 

5 . Subroutine  FIELD 

This  routine  extracts  the  appropriate  bits  of  the  Hollerith 
words  in  DATA  for  each  case  attribute  used  by  the  main  program,  PCP. 

It  then  converts  some  of  these  to  integer  numbers  and  leaves  others 
in  their  original  alphanumeric  form.  Below  is  the  list  of  elements 
extracted  for  a case.  See  Report  10436,  for  the  program  listing  FIELD. 
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15 

12 

II 

14 

18 

12 

12 

12 

12 

12 

12 

A1 

13 

II 

12 

15 

16 

14 


NAME 

DESCRIPTION 

A3A 

MU  controller  code 

A3 

MU  case  no. 

A4A 

Month  notified 

A4B 

Year  notified 

B3 

No.  of  people  on  board  distressed 
(if  B3  > 4095,  B3  = 4095) 

B5 

Value  of  distressed  unit  (if  B5  > 
B5  = 130001) . 

B6 

Air  temperature  (if  unknown,  B6  = 

B8 

Wind  force  (if  unknown,  B8  = 1) 

B9 

Sea  swell  (if  unknown,  B9  = 1) 

BIO 

Visibility  (if  unknown,  B10  = 99) 

B12A 

Type  of  distressed  unit 

B13 

Length  of  distressed  unit 

B16A 

If  B16A  = S,  case  occurred  on  land 

B16 

Distance  offshore 

B17B 

Severity 

Entry  Point  FIELD2  [1=1 , . . . ,CMAX) 

A1A(I) 

District 

A1B(I) 

OP  FAC 

Cl  (I) 

Date/time  notified 

CIA  (I) 

Time  notified  in  minutes  (tenths) 

C5(I) 

No.  of  resources  assisting 

C7A1CI) 

Degrees  ladutude 

C7A2  (I) 

Minutes  latitude 

unit 

130001, 

99) 
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FORMAT 

NAME 

DESCRIPTION 

12 

C7BXCO 

Degrees  longitude 

12 

C7Bs(I) 

Minutes  longitude 

14 

C9(X) 

Date/time  1st  resource  underway 

11 

LGF 

No.  of  logical  fields  in  case 

Entry  Point 

FIELD3  (1=1, .... ,CMAX;  J=1 , . . . ,NO(I) ) 

12 

KLCI.J) 

Assisting  resource  type 

12 

D3A(I  ,J) 

Date  resource  underway 

12 

D3B(I,J) 

Hour  resource  underway 

12 

D3C  (I , J) 

Minutes  resource  underway 

14 

D3(I,J) 

Time  underway  in  minutes  (tenths) 

12 

D4A(I,J) 

Date  resource  on  scene 

12 

D4B(I,J) 

Hour  resource  on  scene 

12 

D4C(I,J) 

Minute  resource  on  scene 

14 

D4(I,J) 

Time  on  scene  in  minutes  (tenths)  (if 
D4(I,J)  = 0,  D4(I,J)  - D3(I ,J) ) . 

14 

D6(I,J) 

Hours  spent  searching  (tenths) 

14 

D8(I,J) 

Total  time  on  case  (tenths  of  hours) 

12 

D9(I,J) 

No.  of  sorties 

12 

Dll (I ,J) 

Assistance  to  personnel 

12 

D12(I ,J) 

Primary  assistance  to  property 

12 

D13(I,J) 

Secondary  assistance  to  property 
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B.  User's  Guide  for  PCP 

Section  A of  Table  1 describes  in  detail  the  layout  of  the  data 
deck  for  PCP.  It  can  be  assumed  that  all  data  must  be  punched  in 
the  designated  card  columns,  right  adjusted. 

Note  that  the  location  data  for  the  current  district  and  the 
OP FAC  location  data  must  be  changed  for  each  district  in  order  to 
exercise  that  district. 

Section  B of  Table  1 gives  the  exact  input  layout  for  the  first 
district.  Note  that  the  card  columns  1-80  are  indicated  in  the  table. 

Figure  5,  the  deck  setup  for  a PCP  run  is  presented  as  a guide 
for  the  NBS  system  UNIVAC  1108. 


Section  A:  Table  1 

Input  Data  to  PCP  (Cards) 

Location  data  for  current  district  being  exercised: 


CC 

FORMAT 

NAME 

Description 

1-  6 

16 

DIST 

current  district  being  exercised 

7-12 

16 

LAT 

latitude  and  longitude  of  district 

13-18 

16 

LONG 

origin  (minutes) 

19-25 

F7.3 

ZBETA 

longitude  correction  factor 

26-31 

16 

XPT 

X,Y  distance  from  origin  for 

32-37 

16 

YPT 

cases  at  OPFAC  with  no  location  data 

38-43 

16 

XLMT 

Upper  limit  of  X,Y  for  current 

44-49 

16 

YLMT 

district 

50-55 

16 

XLOW 

Lower  limit  of  X,Y  for  current 

56-61 

16 

YLOW 

district 

62-67 

16 

XDELTA 

X,Y  distance  from  OPFAC  for  cases 

68-73 

16 

YDELTA 

with  no  location  data,  but  a 

given  OPFAC . 

(All  distances  in  nautical  miles.) 


Table  1 
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TABLE  1 Continued 


Search  speeds 


cc 

FORMAT 

NAME 

DESCRIPTION 

1-80 

2014 

SOA3CI) 

Speeds  of  advance  for  each  type 

1=1,. ..,70 

of  resource  (Dl)  (knots) 

Needs 

Matrix 

CC 

FORMAT 

NAME 

DESCRIPTION 

1-80 

2014 

INEED(I) 

Explicit  need  for  each  type  of  ass is 

1=1, ...,100 

tance  rendered  (Dll,  D12,  § D13) 

District  centroids 

and  probabilities 

(followed  by  EOF) 

CC 

FORMAT 

NAME 

DESCRIPTION 

1-2 

12 

I 

District 

9-10 

12 

DLAT  ( I )") 

Degrees  5 minutes  latitude  of  dis- 

11-12 

12 

MLAT(I)J 

trict  centroid 

18-20 

13 

DLON(I)l 

Degrees  § minutes  longitude  of  dis- 

21-22 

12 

MLON(I)J 

trict  centroid 

29-33 

F5.3 

PROB(I) 

Probability  of  a C-130  case  occurring 
in  district  I,  accumulative. 

Opfac 

location  data 

(followed  by  EOF) 

CC 

FORMAT 

NAME 

DESCRIPTION 

7-11 

15 

LIST (NS) 

OPFAC 

16-21 

16 

XSTA(NS)1 

X,Y  location  of  OPFAC  from  the 

22-27 

16 

YSTA(NS)  J 

district  origin 

NS=1,...,  no.  of  OPFACs  with  location  data 

Random  number  seed 

CC 

FORMAT 

NAME 

DESCRIPTION 

1-12 

012 

B(l)| 

Octal  seed  used  as  input  to  random 

14-25 

012 

B(2)  \ 

no.  generator. 
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C.  Example  PGP  Output  and  Interpretation  of  Output 

Table  29  Output  Data  from  PCP,  describes  each  parameter, 
its  format  and  PCP  name,  as  it  appears  on  the  SARSIM  CASE  FILE 
Figure  2 page  48  is  an  Example  of  a Case  Description  as  pro- 
cessed in  PCP,  with  the  same  format,  sequence  and  description  as 
appears  in  Table  2 . 
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WORD 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 


FORMAT 

12 

13 

15 

13 

16 

14 
18 
12 
12 
12 
12 
12 
12 

11 

12 

II 

11 

12 


Table  2 

Output  Data  From  POP  (Tape) 


NAME 

A1A(1) 

STANO 

A3 

A4 

MINC1 

B3 

B5 

B6 

B8 

B9 

BIO 

B12A 

B13 

NSEV 

NN 


DESCRIPTION 

District 

Station  number 

Case  Number 

Month/Year  Notified 

Date/time  notified 

(minimum  Cl  for  M(J  case) 

No.  of  people  on  board  distressed  unit 
Value  of  distressed  unit 
Air  temperature 
Wind  force 
Sea  swell 
Visibility 

Type  of  distressed  unit 
Length  of  distressed  unit 
Severity 

No.  of  explicit  needs 


(NN=N  if  M=0 ; NN=N+1  if  M>0) 


N 

M 

GAMMA 


No.  of  needs  except  search  and  tow 
No.  of  resources  that  towed  or 
escorted 

Degree  of  non-parallelism  for 
multi-resource  case 
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TABLE  II  Continued 


19 

14 

IE 

Time  on  scene  C tenths  of  hours) 

20 

11 

SI 

No.  of  parallel  long  search,  res. 

21 

12 

S2 

Type  of  short  search  (code) 

22 

14 

TSM 

Total  search  miles  (whole  time) 

23 

15 

MILES 

Distance  offshore  (tenths  of  mi.) 

24 

15 

x ) 

Location  of  case  from 

25 

15 

Y 1 

district  origin 

26 , 

29,.,.. 

12 

NEED  (I) 

Explicit  need  (1=1.. NN) 

27, 

30,.  . . . 

14 

TOS(I) 

Time  on  scene  (1=1.. NN) 

28, 

31,.... 

12 

DELTA  (I) 

Delay  of  next  resourse  relative 

to  previous  resource  (1=1.. NN) 

Note: 

The  above 

data  describes  one  case. 

Each  case  is  written  on 

the  output  tape  in  an  unformatted  binary  write.  This  may  take  one  or 
two  write  statements  in  PCP  depending  cn  the  size  of  the  case.  The 
1st  write  dumps  28  words  on  tape;  the  2nd  write  dumps  the  remaining 
words  needed  to  describe  the  case. 
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1 

jj 

OF  * 

Case  D&aoniofi} 

As 

TRocESSEb 

TCP 

1 

I 1 

5 

20003 

127  52330 

32 

130001 

63 

36 

10 

8 

12 

201 

! 

j I 

1 

0 

0 423 

0 

0 

0 

9 

2 

2 

i 3 

423 

0 

• 

l i 

5 

20001 

78162325 

0 

0 

72 

4 

3 

12 

70 

0 

</ 

i i 

1 

0 

Q 0 

0 

0 

0 

9 

2 

2 

ka 

0 

0 

1 

i 

6 

20001 

116250140 

1 

0 

63 

18 

5 

3 

70 

0 

a 

l 

3 1 

1 

0 

0 810 

0 

0 

0 

9 

2 

2 

46 

n 

010 

0 

j 1 
1 1 

6 

20002 

116261955 

1 

0 

51 

15 

4 

0 

70 

0 

V 

1 

0 

0 484 

0 

0 

0 

9 

2 

2 

til  6 

J 

S 1 

484 

0 

6 

20003 

27210554 

1 

30000 

33 

37 

7 

10* 

60 

0 

T 

S l 

0 

1 

0 0 

1 

0 

9 

9 

2 

2 

.17 

0 

0 

4 

’ 1 

l 

6 

20001 

97290840 

0 

0 

51  * 

20 

6 

10 

70 

0 

H 

1 

1 

0 

0 0 

0 

0 

0 

9 

“7 

66 

31  a 

0 

0 

' 

ii 

6 

20002 

127281658 

0 

0 

66 

1 

0 

10 

70 

0 

r** 

J 

3 1 

1 

0 

0 0 

0 

0 

0 

9 

2 

2 

U8 

4 

0 

0 

i i 

6 

20003 

18  51100 

0 

0 

69 

20 

5 

10 

70 

0 

f 

j i 

1 

0 

0 524 

0 

0 

0 

9 

2 

2 

as 

,i 

524 

0 

J . 
) J- 

6 

20001 

108130830 

0 

0 

49 

10 

8 

10 

70 

0 

-1  1 

1 

0 

0 28 

0 

0 

0 

9 

2 

2 

43 

j 

28 

0 

j i 

7 

20001 

86231600 

1 

0 

57 

18 

3 

10 

70 

0 

c 

i 

1 

0 

0 0 

0 

0 

e 0 

9 

2 

2 

jia 

0 

0 

1 i 

7 

20002 

17  91212 

85 

130001 

99 

25 

1 

99 

11 

201 

f 

a»  l 

1 

0 

0 65 

1 

0 

69 

9 

2 

2 

14 

65 

0 

L i 

7 

20003 

37191601 

1 

75000 

27 

16 

4 

15 

11 

0 

H 

? t 

In 

j 

1 

0 

0 26 

0 

0 

0 

3590 

2 

2 

' 

26 

0 

1 1 

7 

20005 

47132115 ‘ 

0 

0 

99 

1 

1 

99 

70 

0 

f 

j 1 

1 

0 

0 0 

0 

0 

0 

9 

2 

2 

318 

C 

0 

1 1 

7 

20001 

87101210 

0 

V 

0 

52 

20 

5 

99 

70 

0 

V 

j 1 

1 

0 

0 102 

0 

0 

* 0 

9 

2 

2 

lie 

i 

102 

0 

? l 

7 

20003 

18151135 

0 

0 

65 

12 

7 

10 

70 

0 

l 

i i 

1 

0 

0 20 

0 

0 

0 

BO 

2 

2 

Ip  Z 
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IV.  DEMGEN  and  HIST 


A.  Introduction  and  Description  of  the  Programs 

The  final  step  of  the  preprocessing  is  the  ordering  of  the  case 
histories  for  input  to  the  operational  simulator  OPSIM.  There  are  two 
alternative  programs  to  do  this,  DEMGEN  and  HIST.  HIST  creates  a chrono- 
logical file  of  cases,  and  DEMGEN  orders  the  cases  randomly  according  to 
the  user's  direction.  These  final  programs  take  a Coast  Guard  historical 
case  load  tape  containing  the  original  assistance  reports  and  converts  it 
to  an  exogenous  event  tape,  the  CASE  TAPE,  suitable  to  be  used  as  input 
to  OPSIM. 

The  tape  input  to  DEMGEN  contains  the  SARSIM  CASE  FILE,  single  unit 
case  histories  arranged  chronologically  by  OPFAC.  The  OPFACs  are  ordered 
numerically  by  their  identification  numbers.  Multi-unit  cases  next  appear 
on  the  tape  chronologically  arranged.  C-130  cases  appear  last  on  the 
tape,  also  arranged  chronologically. 

Each  case  regardless  of  type  contains  a fixed  length  record  of  28 
FORTRAN  words.  This  record  may  be  followed  by  a variable  length  record 
containing  3(N-1)  FORTRAN  words,  where  N,  number  of  needs,  is  defined  by 
the  15th  word  of  the  first  block.  In  sum,  all  cases  contain  at  minimum, 
a fixed  length  record,  which  may  or  may  not  be  followed  by  a variable 
length  record. 

1.  Program  DEMGEN.  Program  DEMGEN  can  generate  the  CASE  TAPE  in  a 
number  of  ways,  depending  upong  the  specific  arrangement  of  the  options 
available  to  the  user. 

In  general,  DEMGEN  reads  a case  from  the  SARSIM  CASE  FILE  and  determines 

from  the  cases'  historical  occurrence  time, into  which  one  of  8 possible 

"boxes”,  located  on  drum,  the  case  is  to  be  inserted.  ("Boxes"  are 

explained  further  on  in  this  section) . The  one  exception  to  this  process 

occurs  when  one  wishes  to  generate  a CASE  TAPE  according 
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to  more  specific  and  selective  criteria  than  a mere  random  duplication 
or  increase  of  the  SARSIM  CASE  TAPE.  In  'this  event  cases  taken  from  the 
PCP  are  put  into  "auxiliary’  bins"  before  being  inserted  into  one  of 
the  8 possible  boxes.  Mien  all  the  cases  appearing  on  the  SARSIM 
Case  File  have  been  thus  inserted  into  boxes,  DEMGEN  takes  the  given 
[input)  starting  time  (day,  hour  and  minute)  and  selects  the  box 
corresponding  to  that  time.  It  then  randomly  selects  a case  from  the  box 
and  assigns  the  given  starting  time  to  the  case,  as  the  case  occurrence 
time. 

DEMGEN  rearranges  the  order  of  the  words  describing  the  case,  does 
some  minor  calculations , and  writes  the  case  onto  the  CASE  TAPE.  The 
arrival  time  of  the  next  case  is  then  calculated,  the  appropriate 
box  is  chosen  and  another  case  is  randomly  picked.  (The  occurrence 
time  of  this  case  is  taken  to  be  the  arrival  time  mentioned  above.) 

The  same  word  order  rearrangements  and  minor  calculations  are  made 
and  this  case  is  then  written  onto  the  output  CASE  TAPE.  This  process 
is  repeated  until  either  a given  elapsed  time  period  has  been  exceeded 
by  the  last  calculated  arrival  time,  or  a given  number  of  cases  have 
been  written  onto  the  tape. 

The  input  card  deck  consists  of  the  following  cards. 

(a)  An  option  card  which  governs  program  flow 

(b)  Optionally,  the  following  cards: 

(1)  PERCENT  card  which  is  the  percentage  by  which  the  historical 
case  load  is  changed. 

(2)  MAX  CASE  card  which  tells  the  program  how  may  cases  to 
generate  onto  the  output  tape. 

(3)  A NODRUM  card  if  it  is  desired  to  increase  the  historical 


-37- 


caseload  selectively;  e.g.,  the  number  of  severe  cases  by  x percent 
and  the  number  of  cases  involving  fires  by  y percent. 

(c)  The  CALENDAR  cards  which  define  the  time  interval  over  which 
the  historical  cases  occurred.  The  calendar  is  used,  in  part,  for 
the  calculation  of  the  historical  hourly  arrival  rates  for  each  of 
the  eight  boxes . Once  a given  input  tape  has  been  used  and  the 
historical  arrival  rates  calculated,  all  subsequent  runs  using  the 
same  input  tape  can  read  the  arrival  rates  into  core. 

(d)  The  SCENARIO  card  which  is  based  upon  the  Fiscal  Year  1969 
(day  1 is  considered  to  be  July  1,  1968  and  subsequent  dates  are 
assigned  a day  number  relative  to  day  1)  and  which  gives  the  start 

and  end,  day,  hour  and  minute  of  the  first,  and  last  case  respectively, 
that  will  be  put  onto  the  CASE  TAPE. 

(e)  The  PEAK  PERIOD  SCENARIO  card  which  defines  the  number  of 
peak  periods  and  the  start  and  end  day  of  each  peak  period  with  re- 
spect to  July  1,  1968. 

(f)  The  HOLIDAY  SCENARIO  card;  which  contains  the  number  of 
holidays  in  the  scenario  and  the  day  number  (with  respect  to  July  1, 

1968)  of  each  holiday. 

(g)  The  Box  I. D. card  which  gives  the  box  number  (1-8)  for  any 
combination  of  peak  or  non-peak  period,  weekday  or  weekend,  day  or 
night,  the  box  number  identified  with  that  combination. 

(h)  SEED  card  containing  the  double  precision  random  number  "seed". 

(i)  NSCEN  card  which  specifies  the  number  of  versions  to  be  generated. 
A more  detailed  description  of  these  cards  is  given  in  Part  C.  of  this 
section,  "User's  Guide  for  DEMGEN" . Before  going  into  a more  detailed 
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discussion  of  the  program  logic,  a description  should  be  made  of  the 
arrangement  of  the  8 ’ 1)0X65"  on  drum.  When  a case  is  read  from  the 
input  tape  into  core,  the  occurrence  time  of  that  case  is  used  to 
determine  into  which  one  of  the  eight  boxes  the  case  should  be  inserted. 
The  box  numbers  are  given  below: 

Box  1 = Peak  Period  Weekdays  Days 

Box  2 = Peak  Period  Weekdays  Nights 

Box  3 = Peak  Period  Weekends  Days 

Box  4 - Peak  Period  Weekends  Nights 

Box  5 = Non-Peak  Weekdays  Days 

Box  6 = Non -Peak  Weekdays  Nights 

Box  7 = Non -Peak  Weekends  Days 

Box  8 = Non-Peak  Weekends  Nights 

Here  weekends  are  taken  to  mean  both  weekends  and  holidays.  Days  are 
defined  to  be  the  hours  0800-2000  inclusive.  The  boxes  are  nothing 
more  than  locations  on  the  drum.  Each  box  is  a drum  unit  with  a 
starting  address  defined  by  a user  written  1/0  table.  For  the  3 year 
District  1 historical  tape  containing  about  12,000  cases  the 
following  table  shows  the  number  of  FORTRAN  words  allocated  to  each 
box,  and  the  number  of  cases  that  actually  were  put  into  each  box 
in  a specific  run. 


: No. 

Logical  Unit  No. 

Capacity  in  Words 

#Cases 

Defined  Dy  I/O  Table 

1 

10 

100,000 

2060 

2 

11 

40,000 

744 

3 

12 

70,000 

2037 
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4 


13 


504 


30,000 


5 


14 


140,000 


3231 


6 


15 


40,000 


969 


7 


16 


50,000 


1431 


8 


17 


30,000 


399 


For  runs  using  historical  caseloads  from  other  districts,  or  where  the  case- 
load becomes  greater  than  12,000  the  user  should  be  aware  of  the  possibility 
that  the  boxes,  as  defined  presently,  may  have  to  be  redefined  so  that  each  box 
is  large  enough  to  hold  all  the  cases  assigned  to  it.  The  current  drum  alloca- 
tion (8  boxes  and  4 auxiliary  boxes)  reserves  580,000  of  the  688,127  storage 
locations  available  to  the  user. 

The  numbered  paragraphs  below  are  a more  detailed  description  of  the  computer 
logic.  The  reader  should  consult  the  DEMGEN  flowcharts  to  better  comprehend 

the  following  steps. 

(1)  The  input  deck  is  read  into  core  storage 

(2)  Optionally,  the  number  of  days  in  each  box  is  calculated 

(XSWIT  (3)  1)  or,  the  historical  rates  of  occurrence  (X)  are  read  into 

core,  (XSWIT  (3)  =1). 

(3)  Counters  are  initialized.  The  logical  drum  units  (the  boxes) 
are  set  to  their  origin  .and  the  address  pointer  (IADR  (I,J)  contains 
the  address  of  the  ,1th  case  in  the  Ith  box)  of  the  first  case  in 

every  box  is  set  to  1. 

(4)  The  first  28  words  of  a case  are  read  into  core  from  the 
input  tape.  If  the  value  of  the  15th  word  is  0 or  1 then  the 
case  has  been  completely  described.  If  the  value  of  the  15th  word 
is  n,  n > 1,  then  a block  of  (n-1)  *3  additional  words  is  read  into 
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core.  If  it  is  desired  to  create  an  output  tape  which  has  the  same  case 
distribution  as  the  input  tape,  (ISWIT  (1)  =1)  , subroutine  BOX  is  called 
immediately.  However,  if  one  desires  to  create  a case  tape  according 
to  more  selective  criteria  (ISWIT  (1)  = 3)  a call  is  made  to  the 
user  modified  subroutine  SELECT,  before  subroutine  BOX  is  called. 

(5)  Subroutine  BOX  determines  from  the  historical  occurrence  time 
of  the  case,  the  box  number  of  the  box  into  which  the  case  is  to  be 
inserted  plus  the  hour  of  the  day  the  case  occurred.  The  arrays: 

(a)  NC(I,J)  = Number  of  cases  in  Box  I,  occurrence  hour  J, 

(b)  CBQX(I)  = Number  of  cases  in  Box  I , and 

(c)  IADR(I,K)  = starting  address  of  the  Kth  case  in  the  Ith  box 
are  then  updated . 

(6)  The  case  is  then  written  onto  the  logical  drum  unit  corresponding 
to  the  box  determined  in  (5),  in  a single  block,  whose  length  is  the 
total  number  of  words  used  to  describe  the  case  [28  + (n-1)  *3,  n > 0]. 
The  next  case  on  the  input  tape  is  accessed  and  steps  (4) -(6)  are 
re-executed  until  a 99  is  encountered  in  the  first  word  of  a 28  word 
record. 

(7)  At  this  point  every  case  on  the  input  tape  has  been  processed 
and  the  8 boxes  filled.  The  exception  to  this  occurs  when  one  is 

to  increase  the  historical  case  load  by  selective  criteria  (ISWIT (1) 

= 3) . When  this  occurs , cases  meeting  the  selective 
criteria  have  been  stored  in  auxiliary  bins  and  they  must  be  randomly 
selected  from  these  bins  and  inserted  into  the  boxes,  before  the  boxes 
can  be  considered  full.  Subroutine  ADD  performs  this  function. 
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(8)  The  historical  mean  hourly  arrival  rate  is  calculated  for 


each  of  the  8 boxes  according  to: 

A i j s NChj/DAYS^  where  NCh  . is  the  number  of  cases  in  the  Ith 
box,  Jth  hour, 

DAYS^  is  the  number  of  days  occurring  in  the 
time  period  represented  by  the  Ith  box  over  the 
time  period  covered  by  the  input  tape.  Note 
that  the  's  can  optionally  be  read  into  core, 
in  which  case  this  calculation  is  by-passed. 

(9)  If  the  input  caseload  is  to  be  changed  by  x percent 
then  the  A^  's  are  multiplied  by  the  factor  (1  + x/100). 

(10)  At  this  point  DEMGEN  is  ready  to  select  cases  from  boxes , 
give  them  new  occurrence  times,  and  write  them  onto  the  output  tape. 
First, a 4 word  code  record  is  written  on  the  output  tape.  The  first 
case  written  will  have  an  occurrence  time  equal  to  the  scenario  start 
day  and  hour  which  are  input  data,  chosen  with  respect  to  the  Fiscal 
1969  calendar. 

(11)  Using  the  Fiscal  69  calendar  and  the  case  occurrence  time, 
(day  and  hour)  the  box  number  of  the  case  is  determined. 

(12)  The  box  number  and  case  occurrence  hour  determine  X.  .. 

ij 

(13)  The  INUMth  case  in  the  box  is  selected  to  be  written  onto 
the  output  tape  according  to: 

XNUM  - CBOX  (I BOX)  *RN  + 1 


where  I BOX  = box  number 

RN  = a random  number  (0  < RN  < 1) , and 
CBOX  (IBOX)  = number  of  cases  in  box  IBOX. 
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(14)  The  address  of  the  first  word  of  the  INUMth  case  in  box, 

I BOX  is  known,  as  is  the  address  of  the  INUM  + 1st  case.  Subtracting 
the  smaller  address  from  the  larger  will  give  the  number  of  words  used 
to  describe  the  INUMth  case.  The  drum  is  spaced  foiward  to  the  INUMth 
case  and  it  is  read  into  core.  The  ordering  of  the  words  is  changed 
to  conform  to  the  output  tape  format  and  the  case  is  written  onto  the 
output  tape,  with  a Mimimm  of  2 card  images  used  to  describe  me 
case. 

(15)  At  this  point  the  occurrence  time  of  the  next  case  is 
calculated.  The  procedure  is  based  on  the  assumption  that  case  inter- 
arrival times  are  exponentially  distributed  over  periods  of  time  less 
than  one  hour.  A well  known  procedure  ^ is  then  used  to  calculate 
DELTA,  the  interarrival  time. 

Using  the  cumulative  exponential  distribution, 

Prob  {DELTA  < t}  = 1 - exp[~X^.t] , 
a sample  can  be  found  by  setting 

DELTA  = -Ln  [RN]/X^ . , where 

RN  = a uniformly  distributed  random  number  on  (0,1  _] 

If  the  DELTA  generated  in  this  manner  is  not  larger  than  one  hour, 
one  proceeds  to  paragraph  (16) . 

However,  if  DELTA  is  larger  than  certain  values  it  is  modified  by  one 
of  two  procedures  designed  to  eliminate  excessively  large  interarrival 
values.  The  technique  used  is  an  option  of  the  user  (specified  by  ISWIT(8)); 
both  are  outlined  in  the  following  chart. 

^For  example,  F.  Hillier  and  G.  Lieberman,  Introduction  to  Operations 
Research , Holden-Day,  Inc.,  1967.  See  p.  447. 
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GENERATE  RN 
A - -Ln  [RN]/X± 


t = 0 


I 


t = t+1  hour 

j = j+1 

GENERATE  RN 
A = -Ln[RN]/A^ 


A = A + t 


Proceed  to  paragraph  (16) 
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(16)  DELTA  is  then  added  to  the  current  time  (in  minutes)  which 
is  then  converted  to  the  proper  format  (0000-2400) . 

(17)  If  the  time  exceeds  2400  hours  the  day  is  updated.  The  new 
day  and  time  is  checked  against  the  scenario  end  day  and  hour.  Genera- 
tion of  cases  ceases  either  when  ; 

(a)  scenario  end  time  is  exceeded  (ISWTT(IO)  = 0)  or 

(b)  the  number  of  cases  on  the  CASE  TAPE  exceeds  an  input 
variable  MAXCA  (ISWIT(IO)  = 1). 

Otherwise  the  steps  described  in  paragraphs  (11)  - (16)  are  repeated. 

(18)  After  the  last  case  has  been  written  onto  the  CASE  TAPE 
DEMGEN  is  finished.  All  that  remains  is  to  print  out  composite  statistics 
describing  the  scenario  generated.  The  output  is  presented  first 

using  Coast  Guard  station  numbers  as  headings  and  then  repeated  using 
OPSIM  station  numbers.  Then,  if  ISWIT(9)  = 1,  the  output  tape  is  re- 
wound and  the  first  50  cases  which  were  generated  are  printed  out. 

If  the  user  desires  to  generate  more  than  one  version  of  this 
scenario,  the  program  returns  to  paragraph  (7)  of  this  report  and  re- 
peats the  process  until  NSCEN  versions  have  been  written  onto  the 
output  tape,  with  each  version  separated  by  an  end  of  file  nark. 
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2 . Program  HIST 


The  computer  program  HIST  is  a UNIVAC  FORTRAN  V program  which 
creates  a historical  demand,  tape*  to  be  used  as  input  to  OPSIM.  This 
demand  tape  contains  a subset  of  the  cases  existing  on  the  PCP  output 
tape,  the  SARSIM  CASE  FILE,  chronologically  ordered  by  their  historical 
arrival  times.  The  description  of  HIST  given  below  is  divided  into 
3 sections,  input,  the  program,  and  output. 

A . Input 

The  only  input  to  HIST  is  the  output  tape  of  PCP  (i.e.,  SARSIM 
CASE  FILE).  This  tape  contains  in  order,  single  unit,  multi-unit, 
and  C-130  cases.  Single  unit  cases  are  ordered  by  OPFAC  identification 
number,  and  within  each  OPFAC  by  historical  date  and  time  of  arrival. 
Multi -Unit  and  C-130  cases  are  ordered  by  the  historical  arrival 
time. 

B.  The  Program 

The  purpose  of  HIST  is  to  produce  a demand  tape  containing  the 
above  cases,  or  a subset  of  them,  in  chronological  order.  This  is 
done  in  the  following  steps.  The  reader  is  directed  to  the  flow- 
charts of  HIST,  Report  10435. 

(1)  The  first  single  unit  case  belonging  to  the  first  OPFAC  on 
the  SARSIM  CASE  FILE  is  read  into  core.  The  case  is  then  checked 
against  any  scenario  limits . Instead  of  creating  a CASE  TAPE  con- 
taining all  the  cases  from  the  SARSIM  File  ordered  chronologically, 
it  may  be  desirable  to  order  a specific  subset  of  the  SARSIM  caseload; 
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e.g.,  only  those  cases  occurring  during  the  first  three  months  of 
1969,  or  only  those  cases  occurring  during  the  summer  of  1968.  In 
this  case  a culling  routine  (usually  consisting  of  a few  lines  of 
FORTRAN  code)  can  be  inserted  in  the  program  immediately  after  the 
case  is  read  from  the  PCP  tape.  If  the  case  falls  within  these 
limits,  it  is  accepted  and  is  output  onto  drum  and  the  starting  address 
of  the  case  is  noted.  This  starting  address  will  become  the  pointer 
giving  the  address  of  the  first  case  in  the  bin,  into  which  all  single 
unit  acceptable  cases  belonging  to  the  first  OPFAC  are  inserted.  An 
auxiliary  pointer  is  also  initialized  at  this  time. 

(2)  As  succeeding  cases  belonging  to  this  OPFAC  are  accepted, 
they  are  written  sequentially  onto  drum  and  the  auxiliary  pointer 
is  incremented. 

(3)  As  soon  as  an  acceptable  case  is  encountered  bearing  a new 
OPFAC  I.D.  code  a new  bin  is  created  which  has  a starting  address 
equal  to  the  current  value  of  the  auxiliary  pointer. 

(4)  Succeeding  acceptable  cases  belonging  to  the  new  OPFAC  are 
then  written  onto  drum  sequentially  behind  the  first  case,  the  auxiliary 
pointer  being  updated  each  time  a case  is  inserted  into  the  bin. 

(5)  This  process  continues  until  all  single  unit  cases  have  been 
read  from  the  SARSIM  CASE  FILE, culled,  and  inserted  into  their 
proper  OPFAC  bins . 

(6)  Hie  next  type  of  case,  the  multi-unit  case  is  read  from  the 
SARSIM  CASE  FILE,  and  Is  •inserted,  after  culling,  into  the  multi-unit 
bin;  the  current  value  of  the  auxiliary  pointer  at  the  time  the 
first  acceptable  multi -unit  case  is  encountered  becomes  the  starting 

address  pointer  of  the  multi -unit  bin. 
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(7)  The  procedure  in  paragraph  (6)  is  repeated  for  the  last  type 
of  case  on  the  SARSIM  CASE  FILE,  the  C-130  case, 

(8)  When  the  process  outlined  in  paragraphs  (1)  - (7)  has  been 
completed,  n + 2 bins  mil  have  been  created,  where  n is  the  number 
of  different  OP FAC’s  encountered  in  inserting  the  single  unit  cases 
into  bins.  The  starting  address  of  the  first  case  in  each  bin 

the  cases  will  be  ordered  chronologically. 

(9)  The  first  case  in  each  of  the  n+2  bins  is  then  examined, 
the  case  with  the  earliest  arrival  time  is  chosen,  and, after  modi- 
fication of  the  case  descriptors,  is  written  onto  the  output  Case  Tape. 

The  address  pointer  of  the  bin  from  which  the  case  has  been  picked, 

is  incremented  and  now  points  to  the  first  word  of  the  next  case  in 
that  bin. 

(10)  The  process  described  in  (9)  is  repeated  for  the  new  set 
of  n+2  cases,  each  occupying  the  top  position  in  its  bin.  If  at 

any  time  during  this  process  a bin  becomes  empty,  a bin  switch  is  set 
to  eliminate  the  empty  bin  from  consideration. 

(11)  The  program  stops  when  all  bins  are  empty. 

C.  Output 

The  output  of  HIST  is  the  historical  demand  tape,  the  CASE 
TAPE,  a Fortran  formatted  tape,  containing  variable  length  records. 

The  first  record  is  a 4 word  initializing  record  required  by  OPSIM, 
containing  a 3 in  the  first  word  and  zeros  in  the  remaining  three. 

The  last  record  on  the  demand  tape  also  has  4 words,  containing  a 2 in 
the  first  word,  the  day  number  plus  2 in  the  second  word,  where  day  number 
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is  the  day  number  (relative  to  the  day  number  of  the  first  case)  of 
the  last  case  on  the  demand  tape.  The  third  and  4th  words  are  zero. 

The  format  of  the  variable  length  record  describing  each  case  is 
shown  in  Table  4. 
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Table  3 


DEMGEN  and  HIST  Output  Format  of  a Variable  Length  Record  Describing  a Case 

The  record  is  composed  of  a variable  number  of  card  images , a 
minimum  of  two  card  images  is  needed  to  characterize  a case. 


CARD  1 

FORMAT 

(13,  14,  13,  12,  13,  1115) 

Column  No. 

Case  Descriptor 

1-3 

IDNO=l  (the  same  for  every  case  on  the  tape) 

4-7 

Occurrence  day  of  case 

8-10 

Occurrence  hour  of  case 

11-12 

Occurrence  minute  of  case 

13-15 

District  No. 

16-20 

Station  No. 

21-25 

Case  No. 

26-30 

Number  of  People  on  Board 

31-35 

Air  temperature 

36-40 

Wind 

41-45 

Sea  Height 

46-50 

Visibility 

51-55 

Type  of  Vessel 

56-60 

Length 

61-65 

Severity  of  case 

66-70 

N (Number  of  explicit  needs) 

CARD  2 

FORMAT 

(1015,  110) 

Column  No. 

Case  Descriptor 

1-5 

NNN  (Number  of  Needs  other  than  search 

or  tow) 

6-10 

MVM  (Number  of  resources  that  towed  or  escorted) 

11-15 

GAMMA 

16-20 

SIS 

21-25 

S2S 

26-30 

TSM  (Total  Search  Miles) 

-SO- 


Table  3 Continued 


31-35 

36-40 

41-45 

46-50 

51-60 


Miles  (Distance  offshore) 

X 

Y 

Box  No.  (Box  from  which  case  was  selected.) 
Value  of  Vessel 


CARD  3 and  Succeeding  Cards 

If  N (the  number  of  explicit  needs)  is  greater  than  zero  and  MMM  = 0 
then  N values  of  the  group 

NEED(I)  = the  I.D.  number  of  the  Ith  need 
TOS(I)  = Time  on  scene  for  the  Ith  need 
DELTA(I)  = Delay  of  Ith  need 

are  calculated.  These  are  written  onto  the  Output  Tape  in  card  image 
format,  6 groups  to  a card  with  FORMAT  (6  (12,  F6.4,  F4.2)) 

Card  3 NEED(l)  , T0S(1)  , DELTA(l)  NEED (6)  , T0S(6),  DELTA (6) 

Card  4 NEED(7)  ,TOS(7)  , DELTA(7)  NEED(12)  , T0S(12)  , DELTA(12) 

Card  5 


However,  if  N > 1 and  MMM  (the  Number  of  resources  that  towed  or  searched) 
> 0,  then  the  last  group  (NEED(N) , TOS(N) , DELTA (N) ) must  be  written 
on  tape  as  a separate  card  image  with  FORMAT  (12,  F6.4,  F4.2). 
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B„  Subroutines  Within  IBMGEN 


DEMIEN  uses  four  FORTRAN  subroutines , an  assembly  language 
random  number  generator  and  a user  supplied  I/O  table  which  defines 
the  drum  units  and  capacities  of  the  8 'boxes'*  plus  auxiliary  units 
corresponding  to  auxiliary  bins  if  needed.  Descriptions  of  each 
of  these  routines  are  found  in  the  following  paragraphs. 

1.  FORTRAN-  Subroutines 

(a)  Subroutine  BOX 

Subroutine  BOX  is  called  when  a case  is  accessed  from  the  input 
tape.  Given  the  historical  occurrence  time  of  the  case,  BOX  will 
determine  the  number  of  the  box  into  which  the  case  is  to  be  inserted 
along  with  the  hourly  interval  (1-24)  that  the  case  occurred. 

(b)  Subroutine  SELECT 

SELECT  is  called  whenever  XSWlT(l)  - 3,  It  must  reflect  the 
specific  criteria  by  which  a case  is  selected  and  put  into  an 
auxiliary  bin.  Since  the  selective  criteria  may  change  from  run  to 
run,  the  user  must  supply  the  1 or  2 lines  of  FORTRAN  code  which  does 
the  culling.  Two  examples  of  this  are  given  below p and  in  Figure  5. 

Example  A is  an  example  of  the  coding  of  SELECT  for  a computer 
run  approximating  the  SARSIM  caseload,  but  increasing  the  cases  where 
the  number  of  people  on  board  are  greater  than  3,  by  50%.  A single 
auxiliary  bin  is  defined  by  the  I/O  table. 

Example  B is  an  example  of  the  coding  of  SELECT  for  a computer 
run  where  the  number  of  cases  whose  severity  is  greater  than  4 is  to  be 
increased  by  40%,  and  the  number  of  tow  cases  increased  by  70%,  other- 
wise approximating  the  SARSIM  case  load.  Two  auxiliary  bins  are  defined 


by  the  I/O  table. 
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1 ♦ 
2* 

3 * 

4 * 
& * 
b * 
7* 
a* 
v* 

1 o* 

l l * 
i 2 * 
1 3* 
1 4 * 
1 5* 
1 6 * 
1 7* 

l a * 

1 7 * 

2n* 

2 l * 


SUBROUTINE  select 
COMMoN/SEaU/  n up per 

DIMENSION  IADR(B»M500)  , N A D r ( 4 » 1 0 0 U ) , C P E R ( 2 4 ) , C B 0 X ( 8 > ,NCASES(24) 
1 „ I C A S E < 10  o ) 

COMMON  NODKuM,IADR,NADRtNCASES,CPEF?,CBOX»ICASE 

r 

C THESE  ARE  USE*  FURNISHED  STATEMENTS 
C 

I A U N I T s i a 

I F ( I CASE  ( 6 ) # L T o 4 1 GO  TO  1 
J * 1 

NCASES(J)=NCA5E5(J)M 
KsNCASESI J ) + 1 

NADR(J,K)«NADR(J»K-1)4-N  UPPER 
CALL  N j R A N ( !fl,  1 .NuPPER.ICASEJ  1 ) , L ) 

20  I F { L ♦ I > 25  » 20 

1 RETURN 
2 5 .V  R I T E ( 6 , b 0 7 ) 

507  FORMAT  ( 2 X s " ERROR  IN  NTRAN  TRANSMISSION  IN  SELECT") 

STOP 

END 


1 * 
2* 

3* 
4 * 
5* 
b * 
7* 
0* 
9* 
1 0* 
1 1 * 
1 2* 
1 3* 
1 4 * 
15* 
1 6* 
I 7* 
1 0 • 

1 9* 

2 U * 
21  * 
22* 


SUBROUTINE  select 
COM MON /SEAL/  N UPPER 

DIMENSION  lADRf9»450G)lNADR(4,1000)»CPER(24)»CB0X(e),NCASES(24) 
1 » l CASE  ( 1 00  ) 

COMMON  NOqRuMjIADRjNAOR.NCaSE-S.CPER.CBOX.ICASE 

c 

C THESE  are  USE*  FURNISHED  STATEMENTS 

c 

I D R U M s o 

IF(!CA5E(|H)«GT»35  I D R U M * 1 
IF(ICASE(!7),GT.05  I D R U M « 2 
I F ( I DRUM  , E1*)  • 0 ) GO  TO  l 
I AUN I T s | 7 + I DRUM 

N C A S E S ( IDRUM)eNCASES(  I D R U M ) + l 
K ■ N CASES!  I D R u M ) ♦ l 

NADR(  IDRUM»K)*NADR(  IDRUMjK-1  )+NUPPER 
CALL  NTRAN(lAUNlT,l,NUPPER,lCASE(l)  , L) 

20  IF (L+  1 ) 2 5 » 2 0 
1 RETURN 
25  WRITE(6(507) 

507  FORMAT  !2X»*ERR0R  IN  NTRAN  TRANSMISSION  IN  SELECT") 

STOP 


FIGURE  6 
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[c)  Subroutine  ADD 

ADD  is  called  whenever  ISWXT(l)  = 3.  When  the  input  tape 
containing  the  historical  tape  has  been  read,  non-culled  cases  are 
located  in  the  8 bins.  The  culled  cases  are  located  in  the  auxiliary 
bins,  all  cases  of  a specific  type,  in  the  same  bin.  Subroutine  ADD 
tahes  the  cases  out  of  the  bins  in  a random  fashion  and  inserts  them  into 
their  proper  box,  until  the  desired  percentage  increase  for  that  type 
of  case  has  been  reached. 

(d)  Subroutine  ADT1ME 

AUTIME  is  called  whenever  ISWIT(8)  = 1,  and  whenever  the  inter- 
arrival time  of  the  next  case  exceeds  60  minutes.  It  calculates  the 

interarrival  time,  DELTA,  according  to  the  formulae  given  in  Section  IV.A.l. 

2.  The  I/O  Subroutine. 

Below  is  a listing  of  the  current  I/O  user  supplied  subroutine 
which  defines  the  sizes  of  the  8 boxes  plus  four  auxiliary  bins.  The 
auxiliary  bins  have  a capacity  of  25,000  FORTRAN  words  each.  If  more 
auxiliary  bins  are  to  be  used,  the  I/O  table  must  be  modified  in  accor- 
dance with  the  instructions  In  Table  5. 

3.  The  Random  Number  Generator 

This  subroutine  is  an  assembly  language  routine  written  by  the 
Computer  Services  Division  at  N.B.S.  The  user  must  supply  the  seed. 

Any  odd  octal  number  can  be  used.  The  standard  seed  is  the 
number  described  In  Section  IV. C,  User’s  Guide  for  DEMGEN. 
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IOW 


X<J 


Ai 


YOUR  GWM  FORTRAN  I/C 


■ Cl) 


FORTRAN  users  nay  assemble  their  own  I/O  Table  replacing 
the  standard  on  the  EXEC- I I system.  The  I/O  Table  cor- 
relates logical  reac*  and  write  units  referenced  an  a 
FORTRAN  I/O  statement  (such  as  “READ  INPUT  TAPE  1C , ARRAY" 
with  particular  hardware  devices  in  the  computer  room 
(card  reader  , magnetic  tape  umsarvos  , etc.  j * 


/ 


Example  I/O  Table 


7. 

8' 


n: 


ASM 

10, 

10 

+3 

• 

standard 

PUNCH 

+2 

• 

standard 

PRINT 

+ 1 

« 

standard 

CARD  READER 

-4 

* 

unit  0 = 

CARD  REREAD 

-4-2 

* 

UNIT  I = 

PRINTER 

+ 1 

e 

UNIT  2 = 

CARD  READER 

+0 

• 

UNIT  3 (NOT  USED) 

-J- ? A * 

• 

UNIT  4 = 

tape  a 

END 

(Words  after  a 
p 0 X IL  G Cl  "t  X 0 X G 

ci  s c s g gc, 

q r*i  o x g?  o oy  o 1 1 
Assembler  .• ) 


The  signed  values  in  this  Table  each  stand  Tor  a particular 
hardware  device;  the  1 ines  in  this  Table  each  rebate  to  a 
logical  unit  number  or  a standard  unit.  Therefore,  a 
value  entered  on  a particular  line  connects  a hardware 
device  to  a unit  number. 

In  the  more  detailed  explanation  that  follows,  rexer  zo 
the  example  above. 


^This  section  is  reproduced  from  the  NBS  Computer  Services  Division, 
Systems  Notice  #2,  February  26,  1968. 
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Explanation 


1)  The  values  in  this  table  each  stand  for  a particula 
hardware  device;  i.e., 


■t  1 

always 

points 

to 

the 

card  reader 

+2 

always 

po ints 

to 

the 

printer . 

+3 

always 

points 

to 

the 

card  punch. 

-4 

always 

points 

to 

the 

card  reread 

unit  (FORTRAN  V) . 

+'Af  always  points  to  tape  A, 

A complete  table  of  these  values  is  given  later 
in  this  section. 


2)  The  programmer  who  made  this  I/O  Table  used  logical 
units  0,  1,  2,  and  4 in  his  program;  i.e.,  he  wrote 
I/O  statements  such  as  the  one  below  which  references 
unit  2: 

READ  (2,99)  (KARD(I),  I = 1,14) 

99  FORMAT  (13A6,A2) 

In  the  I/O  Table,  he  must  make  entries  as  follows: 

On  line  NTAS$  enter  the  hardware  device  value  for 

unit  0. 

On  the  line  NTAB$  +1  (i.e.  one  line  after  NTA3$) 
enter  the  hardware  device  value  for  unit  1. 

On  line  NTAB$  +2  (i.e0  two  lines  after  NTAB$)  enter 
the  hardware  device  value  for  unit  2, 
etc . 

In  the  READ  statement  shown  above,  the  programmer  intended 
that  unit  2 refer  to  the  card  reader.  Therefore  on  line 
NTAB$  +2 , he  entered  the  value  +1 , which  points  to  the  card 
reader  (as  explained  in  paragraph  1).  The  other  values  were 
entered  similarly:  Note  that  since  unit  3 was  never  referenced 

in  the  program,  he  entered  +0  on  line  NTAB$  +3»  Any  value 
might  have  been  used  here. 
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3)  This  progr amine  r also  referenced  the  standard  PUNCH , 
PRINT,  and  READ  units;  i.e.  he  wrote  I/O  statements  such 
as  the  one  below  which  references  the  standard  card  punch 

PUNCH  I 0 , ALPHA , E ETA , GANNA 

10  FORMAT  (3F10.5) 

For  the  standard  units, ■ he  must- make  entries  in  the 
preceding  NTAE$ 0 The  values  entered  on  these  lines 
fixed  as  follows: 

The  line  NTAB$  -3  is  reserved  for  the  standard 
punch  unit;  therefore  enter  the  value  +3 
on  this  line. 

The  line  NTAB$  -2  is  reserved  for  the  standard 
print  unit;  therefore  enter  the  value  -2  on 
this  line. 

The  line  NTAB$  -1  is  reserved  for  the  standard 
card  reader;  therefore  enter  the  value  rl 
on  this  line. 


-i_  1H6  S 
o.ir  G 
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Table 


Format  of  the  I/O 


As  stated  before,  the  I/O  Table  is  in  1108  Assembler  forma 
Rather  than  define  this  format,  it  is  simpler  to  give  a 
few  rules  relating  to  the  I/O  Table. 

1) 1  NTAB$*  must  start  in  column  1. 

2)  The  asterisk  (*)  following  NTAS$  is  required; 

it  makes  this  label  accessible  to  EXEC- II. 

3) „  Values  may  start  in  or  beyond  column  2. 


4}  At  least  one  space  must  separate  the  label 
NTAB$*  from  the  value  following  it. 

5)  A period  starts  the  comments  field;  at  least  one- 

space  must  precede  the  period. 

6)  The  last  line  must  contain  the  word  ERD  starting 

in  or  beyond  column  2. 
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Hardware  device  values 


The  values  which 
given  below. 

point  to  part 

icular  hardware  devices  are 

+1 

alway 

s points  to  the  card  reader 

+2 

11 

• or inter 

X 

+3 

ft 

card  punch 

4-4 

n 

console  typewriter 

-4 

11 

card  re— ream  unit  ^ r 0 a TiC-un  V) 

-0 

HI 

card  re-read  unit  (FORTRAN  IV) 

+ ' A’ 

n 

logical  tape  unit  A 

+ 'B' 

« 

$r 

logical  tape  unit  £ 

© 

+ TZ' 

»! 

• 

logical  tape  unit  Z 

-4-  * _ * 

IS 

logical  tape  unit  - 

+ ' ( * 

11 

logical  tape  unit  ( 

+DRUM1 

points  to  a 
- Drum  areas 
in  the  I/O 
under  Drum 

. drum  area,  called  DRUM1* 
require  additional  entries 
Table.  See  the  explanation 
assignment s . 

Tape  and  drum  assignments  require  some  further  explanation. 
See  the  following  page. 


VI 1-5 


-59- 


Tap  e a s s i g nine  n t 

The  following  explains  how  the  I/O  Table  connects  the 
unit  number  in  a FORTRAN  I/O  statement  to  the  final  tape  = 
assignment  made  by  the  operator  at  the  console. 

The  example  I/O  Table  on  page  VIII-1  shows  that  the  pro- 
grammer wishes  unit  4 to  reference  a tape  drive.  He  picks 
logical  label  ’A’  and  enters  it  in  the  I/O  Table  on  line 
NTAB$  + 4.  He  also  punches  an  ASG  control  card  for  his 
job  deck,  containing  logical  label  'A*  and  the  tape  reel 
number  , such  as : 

7 ASG  A = 1103 

8 

When  his  job  deck  is  read  in,  the  ASG  control  card  causes 
the  console  typeout : 

MOUNT  1103 

and  the  operator  mounts  tape  reel  1103  on  an  empty  tape 
drive.  When  the  job  comes  up  for  execution,  the  operator 
types  in  the  physical  number  of  the  tape  drive.  If  it 
is  tape  drive  3,  he  types  in: 

/3  1103 

and  program  execution  continues.  i 
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Drum  assignment 


The  FH432  drums  have  approximately  646,000  cells 
available  to  the  user  for  scratch  area  during  execution 
of  his  program.  Address  bounds  of  the  scratch  area  are 
from: 


1300000  to  3777777  octal 
360448  to  1048575  decimal 

The  user  may  treat  this  area  just  as  he  would  treat 
scratch  tapes. 

The  I/O  Table  may  be  constructed  to  take  advantage  of 
drum  scratch  area.  Use  the  entries  +DRUM1  +DRUM2 

etc.,  in  the  1/0  Table  to  point  to  different  drum  areas'. 
For  each  DRUM  entry,  a set  of  seven  locations  must  be 
set  aside  in  the  NTAB$  table.  The  general  form  of  this 
is 


DRUM  (starting  address) 
(ending  address) 
(starting  address) 
+0 
+0 
+0 
+0 

. ^ 

See  example  on  the  following  page. 
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r-  co 


Example 


Construct  an  input/output  table  which 


(a)  assigns  0 as  the  card  re-read  unit 

(b)  units  1,  2,  3 and  5 as  card  reader 

(c)  unit  4 as  card  punch 
(a)  unit  6 as  printer 

(e)  7 and  8 as  two  different  drum  scratch  areas 

each  300000  cells 

(f)  unit  9 as  tape  A and  unit  10  as  tape  3. 


ASM 


* 


10,  10 

+3 

4-2 

+ 1 
-4 . 


DRUM! 


DO 


DRUM2 


DO 


END 


+ 1 
4-1 
4-3 
4-1 
4-2 

+DRUM1 
4- DRUMS 
+ ’A' 

4'B'  ■ 

4-400001 
4-700000 
4-400001 
4 , 4-0 

4-700001 
+1000000 
+700001 
4 , +0 


(Generates  4 lanes  o-.  zeri^ 


(Generates  *-£  ii nsb  or  z^ro 


(See  note  on  next  page.) 
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Note  : 


The  DO  directive  used  above  generates  identical 
lines  of  coding.  Its  format  is: 

DO  count  , line  of  coding 


represents  a blank  space. 

'DO*  appears  in  the  entry  field.  One  or  more 
spaces  follow. 

’count'  is  the  number  of  times  to  repeat  the 
line  of  coding. 

’count'  is  followed  by  one  space  and  a comma. 

The  next  character  after  the  comma  corresponds 
to  column  1 of  the  line  of  coding  desired. 

Thus , DO  4 , + 0 is  equivalent  to 


+ 0 
+ O 
+ 0 
+ 0 
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standard  input /output  table 

Logical  Unit 
0 

1 

2 

3 

4 

5 

6 

7 through  32 

33 

34 

35  and  36 
37  ana  38 
39  and  40 
41  and  42 
43  and  44 
45  cind  4o 


Assignment 

Reread 

Card  Reader 

Printer 

Cara  Puncn 

Con  sore 

~ -y.  ^2  r>  ~ ^ - 

v^d  i.  d Ac  Lv^ci 

Printer 

T ao e s A a n r ouon  L 

Tape  { 

Tape  - 

Entire  Drum  1300000  to 
3777777  (octal) 

tov/er  ha.:  x-^OOooo 
2537777 

Upper  Half  2540000  to 
3777777 

Lower  Third  1300000  to 
2177777 

Middle  Third  2200000  to 
3077777 

Upper  Third  3100000  to 
3777777 
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C.  User’s  Guide  for  DEMGEN 


This  section  discusses  the  program  options  available  to  the  user, 
the  setting  up  of  the  input  deck  the  format  of  the  output  tape  and 
some  examples  of  typical  runs  showing  the  arrangement  of  the  input 
decV  for  each  run. 

1.  Options  Available  to  the  User 

The  values  of  the  array  ISWITCI)  , which  are  punched  on  the  OPTION 
CARD  and  read  into  memory,  govern  the  logic  in  a given  run  of  the 
program.  These  are: 

WORD  Switch  Setting  OPTION 

ISWIT(l)  1 This  run  will  create  an  output 

case  load  similar  in  char- 
acteristics to  the  input 
caseload. 


ISWIT(l) 


ISWIT(2) 


ISWIT(2) 

ISWIT(3) 


3 This  run  will  create  cases 

similar  to  the  input  tape 
with  x types  of  cases  in- 
creased by  CPER(I)  percent 
(1=1, x) 

0 The  caseload  of  the  scenario 
produced  is  equal  to  that 

of  the  input  tape. 

1 The  caseload  is  increased 
by  PER  percent. 

1 Special  validation  option:  his- 

torical Lambdas  are  read  into  core. 
An  output  tape  will  not  be  creat- 
ed. The  Lambdas  of  the  generated 
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WORD 


Switch  Setting 


OPTION 


cases  will  be  printed  out  along 
with  the  historical  Lambdas. 

ISWIT(3)  0 Normal  operation  of  Program. 

Cases  are  generated  using  historical 
Lambdas  and  an  output  tape  is 
created. 


ISWIT(8) 

0 

The  interarrival  time  of  the  next 
case  is  taken  to  be  < 3/X^ . (60) . 

ISWIT(8) 

1 

The  interarrival  time  is  calculated 
according  to  the  formulae  in  Section 
IV.  A.l . 

ISWIT(9) 

0 

The  cases  on  the  Output  Tape  will 
not  be  printed. 

ISWIT(9) 

1 

The  first  50  cases  of  the  Output  Tape  will 

be  printed. 

ISWIT(IO) 

0 

Cases  will  cease  being  generated 
when  the  scenario  end  day  and  hour 
are  exceeded. 

ISWIT(IO) 

1 

Cases  will  cease  being  generated 
when  a given  number  of  cases , 
MAXCA,  has  been  exceeded. 

Note  that  these  switch  settings  are  mutually  independent  and  no  conflict 
in  settings  can  exist. 

2 . The  Input  Deck 

The  cards  making  up  the  input  deck  are  described  below.  They 
are  arranged  in  the  order  listed. 
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CARD  A 


OPTION  CARD 


Column 

Variable  Name 

Fomat 

Value 

1 

ISWIT(l) 

11 

1 or  3 

2 

ISWIT(2) 

11 

1 or  0 

3 

ISWIT(3) 

11 

1 or  0 

4 

ISWITC4) 

11 

Not  Used  (Blank) 

5 

ISWITC5) 

11 

Not  Used  (Blank) 

6 

ISWITC6) 

11 

Not  Used  (Blank) 

7 

ISWITC7) 

11 

Not  Used  (Blank) 

8 

ISWTT  (8) 

11 

0 or  1 

9 

ISWIT (9) 

11 

0 or  1 

10 

ISWIT(IO) 

11 

0 or  1 
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CARD  B 


PERCENTAGE  INCREASE  CARD  (OPTIONAL) 

This  card  is  used  if  ISWIT(2)  = 1.  Otherwise  it  is  omitted. 

Column  Variable  Name  Format  Value 

1-5  PER  F5.0  The  value  of  the 

percent  increase  of 
the  input  caseload 
desired. 

CARD  C 


MAX  CASES  CARD  (OPTIONAL) 

This  card  is  used  if  ISWIT(IO)  = 1 and  is  otherwise  omitted. 


Column 

Variable  Name 

Format 

Value 

1-5 

MAXCA 

15 

The  number  of  cases 
desired  on  the  output 
tape . 

CARD  D 

AUXILIARY  BIN  CARD(S)  (OPTIONAL) 

This  card  is  used  if  ISWIT(l)  = 3 and  otherwise  omitted.  If  used, 
the  run  will  create  cases  according  to  a more  selective  and  specific 
criteria  than  merely  increasing  the  input  caseload  by  PER  percent!/ 


— In  this  case,  subroutine  SELECT  is  used  and  must  be  changed 
to  reflect  the  specific  selective  criteria  by  which  the  user  desires 
to  accept  a case  and  put  it  in  an  auxiliary  bin.  Note  that  the 
auxiliary  bins  are  drum  logical  units  defined  in  the  I/O  Table 
and  their  limits  must  be  large  enough  to  hold  all  the  cases  assigned  to 
them. 
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Column 

Variable  Name 

Format 

Value 

1-2 

NODRUM 

12 

The  number  of  auxiliary 
bins  into  which  the  cases 
meeting  the  selective 
criteria  are  to  be  inserted. 

3-7 

GPER(l) 

F5.2 

CPER(I)  is  the  percent 

8-12 

CPER(2) 

F5.2 

Increase  desired  for  all 

13-14 

CPER(3) 

F5.2 

Those  cases  in  the  Ith 

• • • • 

• 0 • • 

If 

auxiliary  bin. 

Up  to  twelve  bin  percentages  are  defined  on  the  first  card.  If 
additional  bins  are  used, the  13th  bin  percentage  starts  in  column  1 of 
the  2nd  card. 

CARD  E 

PEAK  PERIOD  CARD(S) 

This  card  contains  the  number  of  peak  periods  in  the  time  interval 

covered  by  the  SARSIM  CASE  FILE  (input  tape)  and  the  start  and  end  week 

in  reference  to  the  input  calendar,  (card  F) 

Column  Variable  Name  Format  Value 

1-2  NPP  12  Number  of  peak  periods  in 

the  time  interval  covered 

by  the  input  tape. 

3-4  IYEARS  12  The  number  of  years 

covered  by  the  SARSIM  CASE 
FILE  relative  to  1966;  1966 
is  assigned  the  value  1. 
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Column 

Variable  Name 

Format 

Value 

5-7 

ISWPP(l) 

13 

ISWPP(I)  is  the  start 

8-10 

IEWPP(l) 

13 

week  of  the  Ith  peak 

11-13 

ISWPP(2) 

13 

period  and  IEWPP(I)  is  the 

14-16 

IEWPP(2) 

13 

end  week  of  the  Ith  peak 

period. 


The  first  card  contains  up  to  10  pairs  of  start  and  end  peak 
periods.  If  additional  peak  periods  are  to  be  defined  a second  card 
is  used  containing  10  peak  periods ; the  l|th  peak  period  start  week 
beginning  in  column  1. 

CARD  F 

HISTORICAL  CALENDAR  CARDS 

These  cards  contain  the  array  CAL(I,J,K)  which  for  a given  month  I, 
day  J and  year  K,  will  equal  the  week  number.  For  example,  CAL  (7,1,1)  = 
1.  (July  1st  1966  is  the  start  week).  For  I a weekday  CAL(I,J,K)  is 
even.  For  I a weekend  or  holiday  CAL(I,J,K)  is  odd. 


Column 

Variable  Name 

Format 

Value 

1-6 

CAL (1,1,1) 

16 

There  are  12  entries  per 

7-12 

CAL(2 ,1 ,1) 

16 

card.  The  first  card  contains 

13-18 

CAL(3 ,1 ,1) 

16 

the  values  of  CAL(I,J,K)  for 

19-24 

CAL  (4 ,1 ,1) 

16 

day  1,  year  1,  and  all  months. 
The  2nd  card  is  for  day  2 and 
year  1 . The  31st  card  for 
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67-72 


CAL(12,1,1) 


16 


day  31  year  1.  A cal- 
endar for  a 3 year  histori- 
cal input  tape  will  contain 
93  cards. 


CARD  G 

SCENARIO  LIMITS  CARD 

This  card  contains  the  start  hour  and  day,  end  hour  and  day  of  the 
scenario,  (relative  to  July  1,  = day  1)  If  a specific  number  of 

cases  are  to  be  generated  (ISWIT(IO)  = 1),  the  end  day  and  hour  of 
the  scenario  cannot  be  determined  in  advance  and  are  left  blank. 


Column 

Variable  Name 

Format 

Value 

1-4 

ISCSH 

14 

start  hour  of  the  scenario 

5-8 

ISCSD 

14 

start  day  (relative  to  July 

l)of  the  scenario 

9-12 

ISCEH 

14 

End  hour  of  scenario 

13-1 

ISCED 

14 

End  day  (relative  to  July 

1 of  the  scenario 


CARD  H 

SCENARIO  PEAK  PERIOD  CARD 

This  card  contains  the  number  of  peak  periods  and  the  start  and 
end  day  of  each  peak  period  of  the  scenario.  The  start  and  end  days 
of  the  peak  periods  are  calculated  relative  to  July  1. 


Column  No. 

Variable  Name 

Format 

Value 

1-2 

NSPP 

12 

Number  of  scenario  peak 

periods 

3-5 

NSPP$(1) 

13 

NSPPS(I)  is  the  start 
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Column 

Variable  Name 

Format 

Value 

6-8 

NSPPE(l) 

13 

day  of  the  Ith  peak  period, 
NSPPE(I)  is  the  end  day  of 

27-29 

NSPPS(IO) 

13 

the  Ith  peak  period,  all 

30-32 

NSPPE(IO) 

13 

taken  relative  to  July  1. 

CARD  I 

HOLIDAY  CARD 

This  card  contains  the  number  of  holidays  in  the  time  interval 
covered  by  the  scenario,  or  if  the  scenario  end  hour  and  day  are 
not  known,  the  number  of  holidays  in  any  time  interval  large  enough 
to  contain  the  occurrence  hour  and  day  of  the  last  case  written  onto 
the  Output  Tape.  It  also  contains  the  day  number,  relative  to  July 


1 

of  each  of  the  holidays 

• 

Column 

Variable  Name 

Format 

Value 

1-2 

NOHOLS 

12 

Number  of  holidays  in  scenario. 

3-5 

IHGL(l) 

13 

IHOL(I)  is  the 

6-8 

IH0L(2) 

13 

day  number  of  the 

9-11 

IHQL(3) 

13 

Ith  holiday 
relative  to  July  1 . 

60-62  IHOL(20)  13 

The  day  numbers  of  20  holidays  can  be  put  onto  the  first  card. 
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Succeeding  cards  if  needed  contain  20  holidays  per  card  each  be- 
ginning in  column  1 and  taking  3 columns  with  an  13  format. 

CARD  J 

BOX  I.D.  TABLE  CARD 

This  card  remains  the  same  for  every  run.  It  contains  the 
values  of  the  array,  ILIST(I,J,K) ; 

1=1,  2;  J=l,2;  K = 1,  2;  where: 

1=1  Peak  Period 

1=2  Non-Peak 

J = 1 Weekday 

J = 2 Weekend  or  Holiday 

K = 1 Day 

K = 2 Night 

ILIST(I,J,K)  for  a particular  configuration  of  I,  J,  K gives 
the  box  number  as  follows : 


Column 

Variable  Name 

Format 

Value  Box  No 

1 

I LI ST (1,1,1) 

11 

1 

2 

ILIST(2 ,1,1) 

11 

5 

3 

ILIST(1,2,1) 

11 

3 

4 

ILIST(2,2,1) 

11 

7 

5 

I LI ST Cl, 1,2) 

11 

2 

6 

ILIST(2,1,2) 

11 

6 

7 

ILIST(1,2,2) 

11 

4 

8 

ILIST(2,2,2) 

11 

8 
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CARD  K 


RANDOM  NUMBER  SEED 


A double  precision  random  number  seed  is  read  into  the  array 


B(l) , B (2) . 

The  seed  is 

read  according  to 

an  octal  format. 

Column 

Variable  Name 

Format 

Value 

1-12 

BCD 

012 

154447730601 

14-25 

B(2) 

012 

255751305264 

CARD  L 

LAMBDA  CARDS  (OPTIONAL) 


If  ISWIT(3)  = 1 the  historical  X — Ts  are  not  calculated,  but 
read  into  core  via  these  cards.  The  array  LAMBDA  (I,J)  contains  the 
value  of  X^.  for  box  number  I (1=1,8)  and  hour  J (J=l,24) 


Column 


Variable  Name 


Format 


I- 5 
6-10 

II- 15 
16-20 
21-25 


56-60 


LAMBDA(1,1) 
LAMBDA(2 ,1) 
LAMBDA  (3,1) 


F5.3 

F5.3 

F5.3 


giving  a total  of  16  cards. 
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3.  Output  Description 

(a)  Output  of  the  program  is  in  two  forms: 

The  printed  output  displays 

1.  the  historical  lambdas; 

2.  the  number  of  observations  Cdays)  in  each  of  the  8 boxes 
over  the  time  period  covered  by  the  scenario;  or  if  ISWIT(IO)  = 1, 
the  time  period  ending  with  the  occurrence  time  of  the  last  generated 
case. 


3.  the  of  the  generated  cases. 

(b)  The  Output  Tape  which  contains  the  generated  cases.  This 
tape  is  a formatted  tape  containing  variable  length  records.  It  starts 
and  ends  with  a 4 word  record.  The  format  and  contents  of  the  record 
describing  the  case  is  shown  in  Table  4, 
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4.  EXAMPLES 


Below  are  two  examples  of  possible  runs.  A printout  of  the  input 
card  deck  for  each  of  the  examples  follows. 

(a)  A run  is  desired  where  the  historical  lambdas  are  to  be 
read  from  data  cards , an  output  tape  is  not  to  be  created  ^scenario 
start  time  and  day  are  July  1,  1968,  hour  1 and  cases  are  to  be  generated 
for  135  days  ending  at  2400  hours  of  the  135th  day. 

The  purpose  of  this  run  is  to  compare  the  historical  against  the 
generated  lambdas . 

(b)  An  output  tape  is  to  be  generated  which  will  contain  2500 
cases,  starting  date  August  10,  1968,  hour  1.  The  input  caseload 
tape  is  to  be  increased  by  20% . 
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D.  Users  Guide  for  HIST 


There  is  no  input  data  deck,  the  only  input  to  HIST  is  the 
SARSIM  case  FILE.  The  only  output  of  HIST  is  the  historical  demand 
tape,  the  CASE  TAPE. 
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V.  SELECT 

ft 

SELECT,  a FORTRAN  program,  was  introduced  as  an  inteimediate 
step  between  the  Preprocessor  and  the  Operational  Simulator  to  alter 
the  occurrence  times  of  the  events  on  the  exogenous  event  tape.  Before 
describing  the  possible  uses  of  SELECT,  an  explanation  of  the  run 
deck  will  be  given. 

The  composition  of  the  run  deck  follows: 


The  first  card  is  the  UNIVAC  1108  run  control  card.  The  next  two 
cards  are  assign  control  cards.  The  exogenous  event  tape  created 
by  DEMGEN  or  HIST  should  be  mounted  on  logical  unit  A.  The  new 
exogeneous  event  tape  to  be  created  by  SELECT  should  be  mounted  on 
logical  unit  B with  a write  ring  in.  Following  the  fourth  card, 
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which  is  a control  card  to  call  the  FORTRAN  compiler,  is  the  source 
language  for  Program  SELECT.  Tills  source  language  is  listed  on  the 
following  pages.  Next  is  the  execute  control  card  succeeded  by  the 
one  required  data  card  to  input  IFIRST,  the  first  day,  and  ILAST,  the 
last  day,  under  the  format  (215).  The  FIN  control  card  indicates  the 
end  of  the  run  stream. 

As  mentioned  above,  SELECT  is  an  intermediate  step  to  alter 
occurrence  times  on  the  exogenous  event  tape;  however,  it  may  or  may 
not  be  necessary  to  run  SELECT  before  running  OPSIM.  The  data  for 
each  event  on  the  exogenous  event  tape  contains  three  numbers  (the  day, 
hour,  and  minute)  representing  the  time  at  which  the  event  occurs. 
Program  OPSIM  considers  the  first  day  to  be  day  zero;  if  the  exogenous 
event  tape  has  been  created  such  that  the  first  day  is  day  zero, 

SELECT  may  not  have  to  be  executed.  If  the  tape  was  created  such 
that  the  first  day  is  day  one,  SELECT  should  be  executed  to  change  the 
times.  For  example,  suppose  that  an  exogenous  event  tape  has  been 
generated  for  the  month  of  July,  and  it  considered  the  first  day  as 
day  one.  Program  SELECT  should  be  run  with  IFIRST  = 1 and  ILAST  = 32. 

A new  exogenous  event  tape,  beginning  at  time  zero  and  continuing  to 
midnight  of  July  31,  would  be  generated.  Note  that  both  DEMGEN  and 
HIST  now  consider  the  first  day  to  be  day  zero. 

Another  possible  use  of  Select  is  to  change  the  time  at  which 
Exogenous  Event  ENBSIM  occurs.  If  the  exogenous  event  tape  was  created 
such  that  the  simulation  ends  at  some  arbitrary  future  time  greater  than 
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the  occurrence  time  of  the  last  case,  SELECT  should  be  executed  to 
set  it  equal  to  midnight  of  the  last  day'  of  the  time  duration  being 
simulated.  For  example,  suppose  that  an  exogenous  event  tape  has 
been  generated  for  the  month  of  July,  and  it  considered  the  first 
day  as  day  zero  but  arbitrarily  set  the  time  for  the  end  of  the 
simulation  to  be  three  days  past  the  occurrence  time  of  the  last  case. 
SELECT  should  be  executed  with  IFIRST  = 0 and  I LAST  = 31.  Again, 
both  DEMGEN  and  HIST  now  correctly  set  the  occurrence  time  of 
Exogenous  Event  ENDSIM  to  be  midnight  of  the  last  day  of  the  time 
duration  being  simulated. 

Lastly,  Program  SELECT  can  be  used  to  select  a portion  of  an 
exogenous  event  tape  for  simulation.  For  example,  suppose  an 
exogenous  event  tape  has  been  created  for  the  months  of  January, 
February  and  March  of  1969.  The  user,  however,  wants  to  simulate 
February  only.  SELECT  should  then  be  executed  with  IFIRST  = 31 
and  ILAST  = 59;  a new  exogenous  event  tape,  beginning  at  time  zero 
and  continuing  to  midnight  of  February  28,  would  be  created. 
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Program  SELECT  Source  Language 


dimension  inisooj ,a«so) » b c so > 

10  read  (5,800, end-1000,  f FIRST , ILAST 
WRITE  (6,800)  I E I RST , ILAST 
gOQ  FORMAT  (215) 

WRITE  (a, 802  ) 

002  FORMAT  ( t 3 0 0 Q » ) 

WRITE  (6,902) 

902  FORM  AT  (»  3 0 0 0 8 ) 

iOQ  ikeep  ■ 0 

READ  ( 7 ,805  > UN(1)»I«1,16) 

80S  FORMAT  ( 13,  I *»,  13,  12,  I 3 » 11  15) 

! DA y ■ I N ( 2 ) 

IF  ( IN  ( 1)  * E 3 • 3 ) SO  TO  100 
IF  ( I N ( 1)  , E Q « 2)  GO  TO  200 

IF  ( I N ( 2 ) « L T , IFIRSt  .OR.  In(2)  ,GE,  J L A S T ) 60  TO  ljO 

I N ( 2 ) ■ IN(2)-IplRST 
I END  * I N ( 2 ) 

IKEEP  * i 

WRITE  (8,805)  ( I N « I)  ,I«1  , U) 

WRITE  (6,905)  ( IN(  I ) , I"1  » U> 

90S  FORMAT  ( I H ,13,18, 13,12,13,1115) 

1 1 0 N® I N ( | 6 ) 

READ  (7,810)  U N ( I ) , I * 1 , 1 1 ) 

8 1 0 FORMAT  (1015,110) 

NNN  a I N ( 1 > 

MMM  b I N ( 2 ) 

IF  ( IKEEP  ,EQ.  0)  GO  TO  130 
WRITE  (8,810)  ( I N ( I)  , I*i  , 1 1 ) 

WRITE  (6,910)  ( I N ( I ) , I ■ 1 , 1 I ) 

910  FORMAT  ( 1H  ,1015,110) 

J 30  IF  (N  * LE  * 0)  GO  TO  500 

IF  (N  « G T » i ) GO  TO  ISO 

IF  (MMM  0 6T , I ) GO  To  ISO 

Read  ( 7 ,a i 5 ) I n ( i ) , a < i j 

815  FORMAT  ( I 2 ® F 6 # 8 ) 

IF  ( IKEEP  »EQ9  o ) SO  TO  500 
WRITE  (8,815)  I N ( 1 ) , A ( 1 ) 

WRITE  56,915)  IN(  1 ) ,A<  1 ) 

9 1 5 FORMAT  ( 1 H , I 2 , F6  • H ) 

GO  TO  500 

J50  IF  (NNN  *EQ,  0)  GO  To  180 
INI  ■ N 

IF  (MMM  ,5Ts  0)  INI  « N-I 

READ  (7,820)  (I  N ( I ) » A ( 1)  » B ( I ) , I * l , I N I ) 

820  FORMAT  |4(J2,F6,1»FH,2)| 

IF  (IKEEP  * EQ  • 0)  GO  TO  170 

WRITE  ( 8,820  ) (INI  I)  ,AU  ) ,B<  I ),!■!,  INI  ) 

WRITE  (6,920)  ? IN(I)  ,A(  I)  ,B(  1 ) , I«1  , INI  ) 

920  FORMAT  ( J H ,6«  I 2 e F 6 , R » F **  . 2 ) ) 

J70  IF  (MMM  « EQ « 0)  GO  To  500 
180  READ  ( 7,825  5 I N < J ) 

825  FORMAT  (12) 

IF  « IKEEP  *E9.  0 ) 60  TO  SOO 
WRITE  ( 8 ,825  ) I N ( 1 ) 

WRITE  (6,925)  I N ( 1 ) 

925  FORMAT  « | H ,12) 

500  IF  (IDAV  .LE.  IlaST)  GO  TO  ioO 
200  ILAST  « ILAST  • I F I RST 
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WRITE  (8,030)  JLAST 
030  FORMAT  ( * 2 * , I M , * 

WRITE  (6,930)  I L AST 
930  FORMAT  ( * 2*  , IH»  » 

ENDFIUE  6 
WRITE  (6,935) 

935  FORMAT  C l HO ) 

GO  TO  10 
1000  STOP 

end 


005*  ) 
005*  ) 
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